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1.  INTRODUCTION 


The  relationship  between  road  traffic  characteristics  and  the 
effective  noise  level  at  a given  receiver  location  is  extremely 
complex.  The  characteristics  of  the  sound  emanating  from  a 
particular  roadway  depend  on  parameters  related  to  the  vehicles 
on  the  road,  such  as  vehicle  type,  density,  speed,  acceleration, 
and  directivity  of  radiation.  The  characteristics  of  this 
sound  as  perceived  at  the  receiver  are  affected  by  topological 
considerations  such  as  location,  orientation,  shape,  and  re- 
flectivity characteristics  of  obstacles;  by  atmospheric  and 
ground  absorption;  by  mi crometeorological  conditions;  by 
the  formation  of  shadow  zones;  by  interference;  and  by  other 
effects  specific  to  the  particular  receiver  location.  Finally, 
the  predicted  annoyance  level  depends  on  the  rating  scheme 
applied  to  fluctuations  in  the  sound  intensity.  Accordingly, 
complex  models  for  noise  sources,  sound  propagation  paths, 
and  human  psychology  are  required  if  accurate  predictions  are 
to  be  obtained. 


Accounting  for  as  many  of  these  factors  as  possible  should  be 
the  ultimate  goal  in  modelling  road  traffic  noise.  Neverthe- 
less, in  order  to  keep  the  problem  tractable  without  unduly 
sacrificing  predictive  accuracy,  the  following  major  limita- 
tions have  been  incorporated  in  the  current  implementation  of 
the  computer  model: 

a.  The  computer  model  handles  only  freely  flowing  road 
traffic.  Stop-and-go  traffic  with  a time-dependent 
number  of  vehicles  per  length  of  roadway  and  effects 
of  acceleration  on  the  radiated  sound  power  are  not 
included.  Also,  effects  of  queuing  or  platooning  of 
vehicles  are  neglected. 
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b.  The  computer  model  assumes  a uniform  atmosphere.  Effects 
of  refraction  and  shadow  zones  are  neglected.  Since  the 
sound  pressure  level  at  a distance  of  more  than  a thousand 
feet  from  a source  can  be  strongly  affected  by  micro- 
meteorological  conditions  of  the  atmosphere,  a considerable 
spread  of  actual  data  must  be  expected  around  the  predicted 
values  at  distances  of  more  than  a thousand  feet. 

c.  The  computer  model  regards  all  the  contributions  to  the 
sound  intensity  at  a receiver  as  incoherent  and  adds  the 
intensities  without  consideration  of  possible  phase  rela- 
tions. This  simplification  can  result  in  major  errors 
for  the  sound  pressure  level  close  to  absorbing  or  re- 
flecting surfaces. 

Secondary  limitations  include: 

a.  Simplified  assumptions  for  the  orientation  and  shape  of 
obstacles  in  the  sound  propagation  path, 

b.  first-order  approximations  for  the  effect  of  ground  cover 
on  the  attenuation  of  sound,  and 

c.  consideration  of  single  reflections  only. 

Despite  these  limitations,  the  computer  program  has  yielded 
predictions  of  A-weighted  noise  levels  that  are  generally 
within  ±3  dB  of  those  measured  in  the  field.  Thus,  the 
approach  used  in  deriving  the  mathematical  models  for  high- 
way noise  appears  to  be  valid,  and  one  might  justifiably 
expect  reasonably  accurate  predictions  over  a range  of  high- 
way and  topographical  conditions. 
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2. 


MODEL  OF  NOISE  FROM  FREELY  FLOWING  ROAD  TRAFFIC 


2.1  Parameters  Describing  Fluctuating  Traffic  Noise 

Due  to  its  random  fluctuation,  the  noise  from  road  traffic  is 
best  described  by  statistical  means.  Commonly  used  for  this 
purpose  are  A-weighted  sound  pressure  levels  L that  are  exceeded 
x percent  of  the  time  ("percentile  levels").  The  level  L,_  that 
is  exceeded  50%  of  the  time  and  the  level  L10  that  is  exceeded 
10%  of  the  time  are  frequently  used  as  single  number  descriptors 
of  fluctuating  road  traffic  noise.  Combinations  of  these  levels 
with  the  level  Lgo  result  in  the  traffic  noise  index 

TNI  = Lgo  +4  (L10-L9Q)  + 30  aB  ( A)  , (1) 


and  in  Robinson’s  original  definition  of  the  noise  pollution 
level  [I] 


L 

e 


+ L -L 

10  90 


(2) 


which  is  for  a Gaussian  distribution  of  sound  pressure  levels 
identical  with 


^NP  ^50  + ^ L 10  ^90  ) + 


<L10“L90> 

57 


(3) 


The  energy  mean  level  L and  the  standard  deviation  ar  of  the 
sound  pressure  level  are  employed  in  the  more  general  definition 
of  the  noise  pollution  level: 

LNP  Le  + 2.56  aL  . 


(4) 


When  Eqs.  (2),  (3),  and  (4)  are  equivalent,  which  is  true  in  a 
first-order  approximation  for  noise  from  highway  traffic,  per- 
centile levels  also  can  be  expressed  in  terms  of  the  energy  mean 
level  L and  the  standard  deviation  aT  of  the  sound  pressure 

0 1j 

level : 


(5) 


L = L + 1.28  aT 

1 0 50  L 


(6) 


L - 1.28  aT 

50  L 


(7) 


The  computer  program  calculates  the  energy  mean  level  Lg  and  an 
approximation  for  the  standard  deviation  of  the  sound  pressure 
level,  which  is  based  on  the  formula 


aL  = 4.35  /ln(l+K2)  dB  , (8) 

where  k2  is  the  variance  of  the  sound  intensity  normalized  by  its 
mean.  Since  the  distribution  of  peak  levels  is  of  particular  in- 
terest, the  computer  program  provides  in  addition  a refined  cal- 
culation of  the  level  L . This  calculation  is  not  based  on  a 

1 o 

Gaussian  distribution  of  levels,  but  it  employs  higher  order 
statistical  parameters  of  the  sound  intensity  for  a more  precise 
calculation  of  peak  level  distributions  [2,5]. 


2.2  Energy  Mean  Level 

To  give  the  user  an  overview  of  the  calculations,  we  derive 
briefly  the  mathematical  formulation  for  the  energy  mean  level. 
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We  first  consider  the  sound  emanating  from  a single  vehicle 
located  at  a point  x on  a road  segment  extending  from  xj  to 
x2 . The  sound  intensity  I at  a receiver  located  a distance 
d from  the  road  (see  Fig.  1)  is  given  as 


I = 


I r 2 
ref  ° 

x2+d2 


10 


-D/10 


(9) 


where  I „ is  the  reference  value  of  the  intensity  at  a short 
distance  rQ  (typically  50  ft)  from  a single  vehicle  on  an 
unobstructed  roadway  site,  and  D is  the  uniform  attenuation 
in  dB  of  sound  from  all  points  of  the  road  segment  considered. 
Equation  (9)  is  based  on  the  assumption  of  hemispherical 
radiation  from  a point  source,  for  which  the  intensity  in 
the  diverging  field  decreases  with  the  inverse  square  of  the 
distance . 

Let  us  now  consider  the  situation  in  which  the  road  segment 
xj5x2  contains  a stationary  flow  of  vehicles  having  identical 
noise  characteristics  and  moving  at  essentially  the  same  speed 
The  mean  vehicle  concentration  X in  vehicles/foot,  is  given 
as 


X = 


( vehicles /mile ) 

v- 5280 


(10) 


where  V is  vehicle  speed  in  miles/hour.  The  mean  sound  inten- 
sity at  the  receiver  due  to  all  the  vehicles  on  the  road  seg- 
ment is  then 


I = I „r $ • 10 
ref  0 


-D/10 


X 


x2+d2 


dx 


(ID 
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We  define  the  angle  at  the  observer  (Fig.  1)  as 


a = t an 


-1  x 


d 5 


so  that  Eq.  (11)  is  evaluated  as 


T T 2 ^ A 1 r\~  D/10 

I = I „r„  -T  Aa  10 

ref  o d 


(12) 


where  Aa  is  the  angle  subtended  at  the  receiver  by  the  road 
end  points  x1  and  x . 

For  T different  groups  of  vehicles,  with  each  group  having  an 

intensity  I „ , and  an  average  vehicle  concentration  A,  , 
r*6 1 j t 

we  obtain  the  mean  intensity 


I 


AalO  D//1° 


T 


l 

t=i 


^t  ^ref , t 


(13) 


Field  experience  has  shown  that  the  reference  level  L „ , , 
defined  as  Lref  t = 10  log  Iref  t dB,  has  approximately  a 
normal  distribution  for  certain  large  groups  of  vehicles 
(such  as  passenger  cars  or  heavy  trucks  travelling  at  roughly 
the  same  speed).  Therefore,  rather  than  considering  a large 
number  T of  small  groups  for  all  kinds  of  different  vehicles 
or  for  slightly  different  speeds,  we  consider  a smaller  number 
S of  broader  groupings.  Accordingly,  we  replace  the  summation 
shown  in  Eq . (13)  as  follows: 


T 


l 

t=l 


10 


L „ ,/10 
ref , t 


1 

2 


'L-L 


ref , s 


ref  ,s 


10L/1°  dL, 


(14) 
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X1  ROAD  SEGMENT  X2 


FIG.l  DEFINITION  OF  PARAMETERS  DESCRIBING  A 
ROAD  SEGMENT  AND  A RECEIVER  LOCATION 
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where  A is  the  average  concentration  of  all  vehicles  of  group 
s 

s,  L „ is  the  mean,  and  a „ is  the  standard  deviation  of 
ref , s rel , s 

a normal  distribution  of  reference  levels.  The  expression  for 
the  mean  intensity  from  traffic  with  S groups  of  vehicles  is 
then 


I = 


—r  Aa 
d 


10 


-D/10 


e 

I 

s = 1 


L 


A 10 

s 


ref  , s 


/10 


l(oref,S/4'35)2 


(15) 


In  general,  multiple  road  segments  of  a complex  highway  con- 
figuration differ  in  the  distance  d from  the  receiver  perpen- 
dicular to  the  road,  in  the  angle  Aa  enclosed  at  the  receiver 
by  two  lines  towards  the  ends  of  the  road  segments,  and  in 
attenuation  D of  sound  from  the  road  segments.  Thus,  the 
mean  intensity  of  sound  from  complex  road  traffic  is  given 
by  the  following  double  summation: 


I = 


£ 


road 

segments 


r ' 


d 


Aa 


£ 


A 


-D/10  ^ref^^ 
10  7 u10 


1 

2 


(oref/^ 


35) 


vehicle 

groups 


(16) 


The  energy  mean  level,  in  dB,  is  then  defined  as 

= 10  log  I 


(17) 


2.3  Cumulants  of  the  Sound  Intensity 

Higher  order  statistical  parameters  of  the  sound  Intensity  which 
describe  the  fluctuations  in  noise  are  the  variance,  the  skewness, 
and  the  excess  of  the  intensity  distribution.  These  parameters 
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are  related  to  the  cumulants  <n  of  the  normalized  sound  intensity 
by  [4] 


variance  = k 

2 

3 l 

skewness  = k3/k2 
excess  = k^/k^ 

As  shown  in  Refs.  2 and  3.  the  cumulants  k of  the  normalized 

n 

sound  intensity  from  freely  flowing  road  traffic  can  be  calculated 
from 


,2 


Kn  “ T 


E 


road 

segments 


'n-1 


E 


,2n-l  , . , 

d vehicle 

groups 


n(Lref-D)/1°  1/2 (naref/4 • 35 > 2 
A10n  rei  e rei 


(18) 


where  I is  the  energy  mean  given  in  Eq.  (16),  the  parameters 
d,  A,  D,  L and  are  the  same  as  those  in  Eq.  (16),  and 


C0  = 02-ai 


Cn  = 2H  Cslnol2  c°s2n  1a2  - sina1  cos2n  1al  + (2n-l)Cn_1]  , (19) 


are  functions  of  the  angles  a and  a2  as  defined  in  Fig.  1,  and 
can  be  calculated  by  the  given  recurrence  relation. 

An  essential  property  of  cumulants  is  their  additivity  for 
statistically  independent  processes.  Use  is  made  of  this  property 
in  the  summation  over  vehicle  groups  and  road  segments.  In  addi- 
tion, background  noise  from  sources  other  than  vehicles  on  the 
roads  under  investigation  can  be  considered. 
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When  the  level  of  fluctuating  background  noise  is  normally 
distributed  about  a mean  Lg  with  the  standard  deviation  Cg,  the 
cumulants  of  the  background  noise  intensity  follow  as 


K 


B ,n 


= 10 


n 

10 


( LB+2 


a2/4.35) 


(20) 


The  mean  value  , must  be  added  to  the  mean  intensity  I from 
the  traffic  noise.  The  cumulants  calculated  in  Eqs.  (18)  and 
(20)  can  then  be  added  after  both  sets  are  normalized  with  the 
total  mean  intensity. 

It  appears  from  Eq.  (18)  that  the  distribution  of  sound  power 
output  from  different  cars  within  a vehicle  group  (and  also  the 
uncertainty  about  it)  as  expressed  by  the  standard  deviation 
°ref  level  has  a maJ or  effect  on  higher  cumulants.  For 

example,  if  an  estimate  of  orep  - 2.5  for  homogeneous  passenger 
car  traffic  is  wrong  by  only  5%,  the  relative  error  of  the  fifth 
cumulant  calculated  from 


Ak 


n 


n 


n/2(n-l) [(aref+Aaref)/4. 35] 2 

e 

n/2 ( n-1 ) ( aref/4 . 35) 2 
e 


1 , 


for  n = 5 amounts  to  more  than  40%.  Consequently,  only  a very 
limited  number  of  cumulants  can  be  reliably  used  for  the  subse- 
quent calculation  of  distribution  functions  of  the  sound  inten- 
sity or  of  the  sound  pressure  level. 


2.4  Percentile  Levels 

V 

The  cumulant  k2  is  used  exclusively  in  a first-order  approximation 
for  the  oercentile  levels  L , L and  for  the  noise  pollution 

90  50 


10 


level.  The  approximation  given  with  Eq.  (8)  follows  from  the 
definition  of  the  variance 


K2  = 


I2  p(I)  dl  , 


(21) 


; o 


with  the  transformation  rule 


I2  p(I)  dl  = 


' o 


2 AL/4  .35  / at  \ at  n 

e p( AL)  d ( AL)  , 


(22) 


where  AL  is  the  sound  pressure  level  relative  to  the  energy  mean 
and  p(AL)  is  the  probability  density  of  the  level. 

For  a normal  distribution  of  levels  with  the  mean  mT  and  the 

1j 


standard  deviation  oT , 

-L» 


AL-m. 


L 


p(  AL)  = -1—  e 
/2tt 


-1/2 


c 


L 


(23) 


and  thus 


K 


2mL/4.35  + 1/2  (t^) 


(24) 


When  is  eliminated  from  Eq.  (24)  using  the  normalization  con- 


dition 


k = 1 = e 

i 


mL/4.35  + 1/2 U 


°L  V 


(25) 


one  obtains  the  result  given  in  Eq.  (18). 
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Higher  order  approximations  involving  the  cumulants  k2  through 
k are  used  for  the  calculation  of  the  percentile  level  L.  . 
(With  slight  modifications  of  the  computer  program,  the  same 
approximations  could  be  used  for  the  calculation  of  percentile 
levels  L , L , or  others.)  The  calculation  is  based  on  the 
assumption  that  the  distribution  of  sound  intensities  at  some 
distance  from  a highway  can  be  well  approximated  by  a Pearson 
Type  III  (or  Gamma)  distribution  [2]: 


p(I?)  = 


1 


3 ' f ( q T ) 


q ' ~1  e-I'/3' 


(26) 


The  cumulants  of  this  distribution  are  [4] 


k ’ = 3 ?n  q'  F(n) 
n ^ 


(27) 


By  moving  the  observer  position  close  in  to  the  highway,  one 
would  expect  that  the  above  distribution  of  noise  still  describes 
the  character  of  the  background  noise.  In  addition,  one  has  to 
consider  a new  distribution  of  noise  from  vehicles  near  to  the 
observer.  The  most  important  portion  of  this  noise  distribution 
is  the  distribution  of  peak  intensities,  which  again  may  be 
approximated  by  a Gamma  distribution: 


p(I") 


1 /I"\q"  1 

iwi  Irq  e 


with  the  cumulants 

Kn  _ gnn  qii  F(n) 


(28) 


(29) 
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The  consideration  of  two  independent  distributions,  each  described 
by  two  parameters,  3’,  q’,  and  3",  qM,  leads  to  a low  parameter 
system  which  allows  us  to  account  for  low  cumulants 


k = k ' + k”  . 

n n n * 


(30) 


of  the  total  intensity 


1=1’  + I” 


(3D 


When  the  total  intensity  is  normalized  by  its  mean  value  the  flow 
parameters  can  be  found  as 


whereby 


6'  = i (1  - /T^rs)  , 


(32) 


1 


g”  = ± (1  + /T=rs)  , 


(33) 


q-  = 2 + 1 + 


1 


/1-rs 


(34) 


" = S-  + 1 1 - 


/l-rs 


s 1 2/2 


)- 


(35) 


r = 6 


k -2k, 

3 2 

K +3k  k 

4 2 3 


(36) 


2k2k4  + 12k2k  3-3k * 


k +3<  k 

4 2 3 


S 


(37) 


Equations  (32)  through  (37)  imply  the  conditions  <3  _>  2k^  and 
rs  < 1,  which  are  commonly  met  for  noise  from  freely  flowing 
road  traffic . 

The  probability  density  of  the  total  normalized  intensity  I 
follows  from  the  convolution  integral  of  the  independent  proba- 
bility densities  p(I’)  and  p(I").  The  result  can  be  written 
with  the  confluent  hypogeometric  function  M [5]: 


p(I) 


q'+q"-l  e-I/3f  n 

^ M[q",q'+q",  I(^-  - g^r 

(3*  )Q  (3")q  r(q'+q")  * ^ 


)]  . 


(38) 


The  cumulative  distribution,  that  allows  for  the  calculation  of 
percentile  levels,  follows  from  Eq.  (38)  by  integration  by  parts. 


I 

p(I)  dl  = 

o 


/ 3 ’ \Q"  v r ( q"+n ) 
Wl  nt0  rTVTnT 


P ( q ’ +q" , 


3 


J)  -e“I/e! 


n-1 

l 

m=0 


(1/3  i ^q'+q’^n-m-l 
f ( q ’ +q"+n-m) 


(39) 


f denotes  the  complete  Gamma  function,  P the  incomplete  Gamma 
function  [£],  and  the  second  term  in  the  brackets  is  used  for 
the  numerical  evaluation  of  incomplete  Gamma  functions  with  the 
additional  argument  n. 
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2.5  Consideration  of  Frequency  Spectra 

In  dealing  with  statistical  distributions  of  sound  pressure  levels, 
it  is  common  practice  that  the  levels  are  A-weighted  levels  at  the 
receiver.  Consequently,  the  cumulants  describing  the  statistical 
distribution  must  be  expressed  in  terms  of  A-weighted  sound  pres- 
sure levels. 

(L  -D)/1° 

In  Eq.  10,  the  level  appears  in  the  expression  10  . For 

rough  calculations,  one  can  assume  that  is  the  A-weighted 

level  at  a reference  distance  rQ  from  a single  vehicle  and  that  D 
is  the  sound  attenuation  in  the  frequency  band  around  500  Hz, 
which  is  most  typical  for  the  attenuation  of  A-weighted  noise  from 
road  traffic.  For  a refined  calculation,  one  must  consider  the 
frequency  dependence  of  the  attenuation  D.  The  computer  program 
allows  for  this  refinement  and  calculates  the  expression 

(Lref-D)/1° 

10  in  octave  bands  centered  at  frequencies  ranging 

from  63  Hz  to  8000  Hz.  A summation  over  the  expressions  in  all 
octave  bands  is  performed  in  order  to  obtain  cumulants  and  sub- 
sequently percentile  levels  that  are  referred  to  the  A-weighted 
noise  at  the  receiver. 

In  the  remainder  of  this  section,  we  describe  the  calculation  of 
the  quantities  employed  in  Eq.  (16)  and  discuss  the  data  used  in 
the  computer  program  for  the  numerical  evaluation  of  the  energy 
mean  level,  the  noise  pollution  level,  and  percentile  levels. 

2.6  Reference  Levels 

The  parameters  L „ and  a r defining  the  normal  distri- 
bution  of  levels  at  a distance  rQ  from  single  vehicles  of 
group  s generally  depend  on  the  velocity  of  the  vehicles 
and  on  the  frequency  band  considered. 
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There  is  no  theoretical  limit  to  the  number  of  vehicle  types 
that  may  be  considered  by  the  program.  However,  the  current 
implementation  provides  data  blocks  for  only  two  types:  cars 

(Type  1)  and  trucks  (Type  2). 


Reference  levels  for  cars  (L  „ , ) are  stored  in  the  computer 

ref , 1 ^ 

program  and  are  shown  in  Pig.  2.  A-weighted  octave  band  levels 
as  well  as  A-weighted  overall  levels  are  given  for  vehicle 
speeds  of  30  and  70  mph.  The  data  are  based  on  measurements 
made  by  N.  Olson  [7].  They  show  the  predominant  influence  of 
tire  noise  with  maximum  levels  at  1000  Hz  for  30  mph  and  at 
2000  Hz  for  70  mph.  An  increase  of  the  A-weighted  overall 
level  by  14  dB  for  an  increase  of  the  velocity  V from  30  to 
70  mph  corresponds  closely  with  the  frequently  observed  speed 
dependence  of  40  log  V.  For  velocities  different  from  the 
stored  data,  the  computer  program  interpolates,  using  the 
formula 


Lref,l(V) 


L „ , (30 
ref  ,1 


mph ) + 


10  log 


_V_ 

30 


10  log 


70 

30 


[Lref,i(70  mPh) 


- Lref,l^30  mPh)J 


(40) 


Below  1000  Hz,  the  octave  band  levels  at  50  mph  are  about  the 
same  as  those  found  by  Galloway  et  al  [5].  Above  1000  Hz, 
Galloway’s  data  are  lower  by  5 to  10  dB. 


In  contrast  to  our  speed  dependence  with  38  log  V for  the  A- 
weighted  overall  level,  Galloway  found  only  30  log  V.  However, 
the  differences  are  only  1.2  dB  at  70  mph  and  1 . 8 dB  at  30  mph. 
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The  standard  deviation  a „ , is  based  on  Olson’s  data.  Although 

ref,l  ° 

there  is  some  dependence  on  frequency  and  velocity,  we  propose  to 

use  the  single  number  a _ = 2 . 5 dB . This  number  is  consistent 

ref,l 

wi'th  findings  in  the  USA  [5]  and  in  England  [5].  Note  that  any 
number  smaller  than  3 dB  increases  the  energy  mean  level  by  less 
than  1 dB . 

Reference  levels  for  trucks  (L  ^ are  plotted  in  Pig.  3 • The 
data  are  based  on  measurements  made  by  Olson  [7]  and  on  recent 
measurements  at  the  New  Jersey  Turnpike.  We  ignore  the  speed 
dependence  observed  in  the  higher  frequency  bands  for  trucks  with 
cross-bar  tires,  and  we  account  for  engine,  fan,  and  exhaust  noise 
exclusively  with  speed-independent  levels. 

The  A-weighted  overall  level  is  higher  by  6 dB  than  data 
reported  by  Galloway  [3].  We  found  out  that  Galloway's  data 
are  referred  to  trucks  with  good  mufflers,  but  that  our  data 
are  more  representative  for  actual  truck  noise.  The  influence 
of  mufflers  on  the  A-level  indicates  that  most  of  the  truck 
noise  is  radiated  from  the  exhaust.  The  effective  source 
height  of  trucks  should  therefore  be  chosen  about  8 ft  higher 
than  the  source  height  of  passenger  cars. 

A value  of  3-5  dB  for  the  standard  deviation  a „ „ is  con- 

ref, 2 

sistent  with  data  from  Galloway  [3],  except  for  steep  gra- 
dients of  the  road  where  considerably  higher  numbers,  up  to 
6.8  dB,  have  been  measured. 
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OCTAVE  BAND  CENTER  FREQUENCY (Hz) 


2 . 7 Attenuati on 


2.7.1  Atmospheric  absorption 

For  outdoor  sound  propagation,  the  attenuation  by  atmospheric 
absorption  at  a temperature  of  68°F  can  be  reasonably  well 
approximated  by  [7(7] 


Da  = 7.4  ■~r~  10"  8 dB 
A cp 


(41) 


where  f = geometric-mean  frequency  of  band,  Hz 
cf)  = relative  humidity  of  the  air,  % 
r = distance  between  source  and  receiver,  m. 


The  computer  program  uses  a relative  humidity  of  55%  and  octave 
bands  with  the  numbers  n = 2,..., 9 for  mid-frequencies  from 
62.5  Hz  to  8000  Hz.  Distances  are  given  in  feet.  Thus, 

Eq . ( 4l ) bee omes 


7-4 


Ln  1000\ 

l2  -&-) 


(l0~8\  r 1 

\ 55  ) 3.28 


10 


- 7 


(42  ) 


where  r’  is  taken  by  the  computer  program  as  the  minimum  dis- 
tance, in  feet,  between  the  receiver  and  a road  segment.  At 
a frequency  of  500  Hz,  where  n = 5,  the  attenuation  is  less 
than  1 dB  up  to  distances  of  10,000  ft.  Since  the  attenuation 
of  traffic  noise  is  typically  characterized  by  the  attenuation 
in  the  500-Hz  octave  band,  atmospheric  attenuation  has  hardly 
any  influence  in  all  practical  cases  of  noise  prediction  within 
a distance  of  a few  thousand  feet  from  a road.  Consequently, 
no  checks  are  made  in  the  computer  program  for  small  differ- 
ences in  attenuation  by  atmospheric  absorption  of  sound  from 
different  locations  on  a road  segment;  it  is  assumed  that  the 
attenuation  from  the  nearest  point  on  the  road  is  representa- 
tive for  the  entire  road  segment. 
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Note  that  the  attenuation  predicted  by  Eq.  (42)  will  be  too  low 
if  the  actual  humidity  of  the  air  is  considerably  less  than  55% 
or  if  the  actual  temperature  is  between  20  and  30°F.  Conversely, 
the  predicted  D will  be  too  high  for  relative  humidities  much 

.H. 

greater  than  55%  \_10~\. 

2.7 .2  Diffraction 

Diffraction  of  sound  is  caused  by  obstacles  in  the  direct  or 
reflected  propagation  paths  from  the  roadway  to  the  receiver. 

Such  obstacles  can  be  artificial  barriers,  earth  berms,  hills, 
buildings,  etc.  For  the  calculation  of  attenuation  by  diffrac- 
tion, we  assume  that  the  obstacle  can  be  modeled  by  a rigid, 
impervious  screen  oriented  perpendicular  to  the  ground  plane 
and  that  sound  is  diffracted  over  the  top  edge  of  the  screen 
exclusively.  We  neglect  the  shape  of  hills  and  the  thickness 
of  barriers  because  of  the  lack  of  available  knowledge.  The 
sound  absorption  and  transmission  properties  of  barriers  are 
not  considered  because  they  play  a minor  role  in  most  practical 
cases.  We  expect  that  the  neglect  of  diffraction  around  the 
ends  of  barriers  will  introduce  no  significant  errors,  and 
it  simplifies  considerably  the  computational  procedures. 
Furthermore,  a diffracting  barrier  is  now  completely  specified 
by  the  coordinates  of  the  two  end  points  of  the  top  line. 

The  attenuation  of  sound  by  barriers  is  determined  primarily 
by  the  difference  6 between  the  path  length  of  the  shortest 
ray  from  the  source  over  the  top  edge  of  the  barrier  to  the 
receiver  and  the  path  length  of  the  direct  ray  from  the  source 
to  the  receiver  in  the  absence  of  the  barrier  (see  Fig.  4). 
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FIG. 4 DEFINITION  OF  THE  PATH  LENGTH  DIFFERENCE  6 FOR 
SOUND  DIFFRACTION  BY  A BARRIER 
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For  large  path  length  differences,  the  attenuation  in  the  acous- 
tical shadow  zone  of  a barrier  is  limited  by  effects  of  refraction 
and  scattering  of  sound  in  the  atmosphere.  Based  on  data  from 
Rath§  [12],  the  computer  program  has  a built-in  maximum  attenua- 
tion of  24  dB. 


The  attenuation  is  not  zero  for  zero  path  length  difference 
(i.e.,  for  a ray  grazing  over  the  barrier).  For  this  situation, 
the  theory  of  Fresnel  diffraction  yields  an  attenuation  of 
about  5 dB.  The  attenuation  becomes  negligible  when  a direct 
sound  ray  travelling  from  the  source  to  the  receiver  passes 
far  over  the  top  edge  of  the  barrier.  To  simplify  computations, 
we  neglect  diffraction  effects  when  the  height  difference  between 
the  direct  ray  and  the  top  of  the  barrier  is  greater  than  20  feet. 

For  height  differences  of  20  feet  or  less,  we  calculate  the 
Fresnel  number 

N = ^ , (43) 


where  6 is  the  path  length  difference  and  A 
corresponding  to  the  center  frequency  of  an 
assume  normal  atmospheric  conditions,  with 
of  1120  ft/sec.  Thus,  for  a center  frequen 
number  is  taken  to  be 


is  the  wavelength 
octave  band.  We 
a speed  of  sound  c 
cy  f,  the  Fresnel 


2f 


(44) 


N 


r> 


6 


We  employ  the  Fresnel  number  and  an  analytic  approximation  to 
measured  data  from  Maekawa  for  the  calculation  of  the  barrier 
attenuation  \_10~] 


j 


(20  log 


/2ttN 


tanh  /2ttN 


+ 5 dB  for  N > -0.2 


DB  ■ 


05) 


■0 


otherwise 


The  formula  Eq . (45)  is  applicable  to  both  positive  and  nega- 

tive values  of  N.  However,  for  the  actual  computation  we 
separate : 


D 


B 


log 

s 

t= 

OJ 

tan/27r  | N | 

log 

/2ttN 

tanh  /2ttN 

dB 

5 dB 


5 dB 


for  N £ -0.2 

for  -0.2  < N £ 0 

for  0 < N £ 12.5 
for  N > 12.5 


(46) 


The  last  line  in  Eq.  (46)  accounts  for  the  above-mentioned 
upper  limit  to  barrier  attenuation. 


As  shown  in  Fig.  55  the  attenuation  of  the  A-weighted  sound 
pressure  level  of  typical  passenger  car  noise  (see  Fig.  2) 
is  almost  identical  with  the  sound  attenuation  in  the  500  Hz 
band.  Hence,  the  only  number  important  for  the  attenuation 
of  road  traffic  noise  is  the  path  length  difference. 
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BARRIER  ATTENUATION  AS  A FUNCTION  OF  THE  PATH  LENGTH  DIFFERENCE 
FOR  A FREQUENCY  OF  500  Hz  AND  FOR  A -WEIGHT ED  LEVEL  OF  TYPICAL 
PASSENGER  CAR  NOISE 


I 


The  path  length  difference  accounts  for  heights  and  distances 
of  a point  source,  a receiver,  and  the  top  edge  of  a barrier. 
Furthermore,  It  allows  us  to  account  for  the  reduced  attenuation 
of  rays  oblique  to  the  top  edge  of  the  barrier  \_12~\. 

For  noise  from  a road  segment  and  for  a barrier  at  oblique 
angle  to  the  road,  the  computer  program  finds  the  path  length 
difference  6^  for  sound  from  the  nearest  point  on  the  road. 

Then,  by  assuming  a monotonic  variation  of  the  path  length 
difference  from  other  points  on  the  road,  the  extreme  ends 
of  the  road  segment  are  considered.  If  the  path  length 
differences  and  62  for  these  end  points  differ  from  6^ 
by  more  than  a number  that  results  in  an  attenuation  differ- 
ence of  about  1 dB,  the  road  segment  between  the  near  point 
N and  the  point  1 or  the  point  2,  respectively,  is  cut  in  half. 
New  path  length  differences  are  calculated  for  the  new  end 
points  of  the  road  segment,  and  the  procedure  of  reducing  the 
length  of  the  road  segment  is  repeated  until  the  attenuation 
by  diffraction  is  approximately  constant. 

The  criterion  used  for  the  acceptance  of  a sufficiently  small 
difference  in  path  length  differences  is  : 

6,+62  / 61+62 

100  l1  + 2 

The  numbers  are  based  on  a frequency  of  500  Hz,  for  which  the 
effect  of  Eq . (47)  on  the  attenuation  is  plotted  in  Fig.  6. 

In  case  of  multiple  diffraction  by  several  barriers  in  parallel, 
we  account  for  the  strongest  diffraction  exclusively.  This  is 
a conservative  procedure  resulting  in  attenuations  that  are 
somewhat  too  small,  but  it  seems  to  be  the  most  reasonable 
way  to  bypass  the  very  complicated  and  not  yet  fully  under- 
stood problem  of  multiple  diffraction. 


< 0.1 


(47) 
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FIG. 6 BARRIER  ATTENUATION  AS  A FUNCTION  OF  THE  FRESNEL  NUMBER  OF  DIFFRACTIONS. 

RESULTS  OF  MAEKAWA'S  MEASUREMENTS,  EQ  . (46  ) WITH 

TOLERANCES  EQ.(47)  AT  500  Hz. 


I 

j 

2.7.3  Ground  absorption 

1 

Ground  attenuation  is  a function  of  the  structure  and  the 
covering  of  the  ground,  both  of  which  influence  its  acoustic 
properties,  and  of  the  heights  of  the  source  and  receiver 
above  the  ground. 

In  the  computer  program,  we  consider  in  a very  simple  approxi- 
mation rectangular  ground  strips,  which  are  defined  by  two 
end  points  of  a center  line  and  by  a width,  and  which  have 
either  a low  cover,  such  as  shrubbery  and  thick  grass,  or  a 
high  cover,  such  as  trees. 

The  height  of  a sound  ray  travelling  from  the  source  to  the 
receiver  over  the  ground  strip  is  checked  only  writh  respect 
to  the  center  line  of  the  strip.  Thus,  it  is  assumed  that 
the  plane  of  the  ground  strip  is  approximately  parallel  to 
a plane  defined  by  a road  segment  and  a receiver.  If  the 
height  of  the  direct  sound  ray  from  the  source  to  the  re- 
ceiver is  more  than  10  ft  above  a ground  strip  with  a low 
cover  or  more  than  30  ft  above  a ground  strip  with  high 
cover,  any  sound  attenuation  due  to  ground  absorption  is 
neglected.  The  heights  of  10  and  30  ft  are  based  on  rough 
estimates  rather  than  on  field  experience  and  might  be 
revised  if  found  necessary. 

In  general,  the  amount  of  ground  attenuation  cannot  be  stated 
in  terms  of  excess  attenuation  per  unit  of  distance.  To  a 
first  approximation,  however,  such  behavior  can  be  assumed 
in  the  range  of  distances  of  200  to  2000  ft  unless  the  total 
attenuation  exceeds  30  dB. 
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No  attempt  has  been  made  to  calculate  accurate  distances  over 
a ground  strip  with  the  computer  program.  Instead,  mean  path 
lengths  r over  ground  strips  are  calculated  with  the  formula 


r 


tt/2 


(48) 


where  w is  the  width  of  the  strip  and  £ the  length  of  the 
center  line. 

The  following  analytical  approximations  to  average  values  of 
measured  data  [20]  are  used  to  calculate  the  attenuation  of 
sound  propagating 

a.  through  shrubbery  and  over  thick  grass 


(see  Fig.  8),  where  r is  given  in  feet  and  f in  Hz. 

For  the.  computer  program,  we  calculate  the  following  attenua- 
tions with  r in  feet  and  with  octave  band  numbers  n,  where 
n = 5 for  the  octave  band  center  frequency  f = 500  Hz: 

a.  for  shrubbery  and  thick  grass 


Dq  = [0.18  log  (f)  - 0.31]  dB 


(49) 


(see  Fig.  7) 


b.  through  tree  zones 


Dq  = 0.01  (f)1/3  dB 


(50) 


! 


(0.0l6n  - 0 . 02 8 ) r dB  if  £ 30  dB 

(jT 
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G 


30  dB 


if  Dq  > 30  dB 


(51) 
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FIG. 7 


FIG. 8 


OCTAVE  BAND  CENTER  FREQUENCY  (Hz) 

ATTENUATION  FOR  SOUND  PROPAGATION  THROUGH  SHRUBBERY  AND  OVER 
THICK  GRASS,  MEASURED  DATA  AND  ANALYTICAL  APPROXIMATION 
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OCTAVE  BAND  CENTER  FREQUENCY  (Hz) 

ATTENUATION  FOR  SOUND  PROPAGATION  IN  TREE  ZONES,  MEASURED 
DATA  AND  ANALYTICAL  APPROXIMATION  FOR  AVERAGE  U.S.A.  FORESTS 


30 


for  tree  zones 


b . 


2n/ 3 
1,310 

30  dB 


r dB 


if  Dq  < 30  dB 
if  Dq  > 30  dB 


(52) 


The  notation  in  the  second  lines  of  Eqs.  (51)  and  (52)  means 
that  the  attenuation  is  limited  to  30  dB. 


Consistent  with  the  assumption  that  the  ground  attenuation  is 
proportional  to  the  mean  path  length  of  the  sound  over  the 
ground  strip,  the  computer  program  accumulates  attenuations 
of  various  ground  strips  in  series.  However,  since  the  path 
length  considered  represents  a statistical  average  for  rays 
propagating  in  all  directions  over  the  strip,  the  path  length 
over  two  equal,  parallel  strips  is  not  just  twice  the  path 
length  over  a single  strip  of  twice  the  width  but  is  generally 
shorter.  Consequently,  the  attenuation  calculated  for  one 
strip  will  be  smaller  than  the  attenuation  calculated  for 
two  strips  in  parallel  each  having  half  the  width. 


The  purpose  of  using  the  statistical  formulation  for  path  length 
given  by  Eq.  (48)  is  to  obtain  reasonable  predictions  for  the 
effects  of  ground  absorption  on  the  average.  There  exist, 
however,  particular  cases  where  the  model  will  not  be  very 
accurate.  For  example,  the  attenuation  of  sound  from  a short 
road  segment  by  long,  narrow  absorbing  strips  is  overestimated, 
whereas  the  attenuation  by  a wide  strip  oriented  perpendicu- 
larly to  the  road  is  underestimated.  To  some  extent,  we 
expect  these  modeling  errors  to  compensate  for  one  another 
in  most  practical  situations.  In  general,  inaccuracies  are 
inherent  to  the  entire  problem  of  ground  absorption. 
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2.8 


Ref 1 ecti on 


The  sound  field  at  a receiver  results  from  contributions  of 
direct  (or  diffracted)  and  reflected  rays.  In  many  practical 
cases  of  sound  propagation  from  a highway,  corrections  applied 
for  reflections  are  small  compared  to  the  inaccuracies  involved 
in  the  prediction  of  ground  attenuation  and  in  uncertainties 
with  acoustical  shadow  zones  owing  to  wind  and  temperature 
gradients  in  the  atmosphere.  Therefore,  the  computer  program 
has  been  designed  to  account  for  reflections  with  a first- 
order  approximation. 

In  the  computer  program,  we  disregard  phase  relations  between 
the  various  contributions  and  consider  incoherent  waves  for 
which  the  total  sound  intensity  is  made  up  by  the  sum  of  the 
intensities  of  the  individual  contributions. 

There  are  always  reflections  from  the  road  surface.  However, 
the  contributions  from  these  reflections  are  already  accounted 
for  in  reference  levels  at  a short  distance  from  individual 
vehicles  on  the  road.  They  don’t  enter  again  into  the  cal- 
culations . 

Reflections  at  the  ground  plane  further  out  from  the  roadway 
are  disregarded  because  they  generally  result  in  a complex 
interference  pattern  with  the  direct  ray.  Consideration  of 
these  effects  is  beyond  the  scope  of  a first-order  approxi- 
mation. 

Reflections  at  any  inclined  plane  result  in  rays  directed 
either  towards  the  ground,  and  thus  being  neglected,  or 
towards  the  sky,  and  thus  not  contributing  to  the  sound 
intensity  at  a normal  receiver  location  close  to  the  ground. 
Therefore,  we  consider  only  reflections  on  planes  that  are 
perpendicular  to  the  ground  plane. 


Within  the  first-order  approximation,  we  further  neglect  the 
actual  frequency-dependent  magnitude  of  reflection  coeffi- 
cients and  distinguish  only  between  reflection  coefficients 
0 and  1 of  reflecting  surfaces. 

In  order  to  determine  whether  a reflective  barrier  is  high 
enough  to  be  effective,  we  consider  the  possibly  reflected 
ray  that  travels  a minimum  distance  from  the  road  segment 
to  the  receiver.  A reflective  plane  perpendicular  to  the 
ground  is  considered  high  enough  if  the  ray  strikes  the 
barrier  at  least  two  feet  below  the  top  edge  of  the  barrier. 

For  reflection  points  within  two  feet  of  the  top  edge, 
diffraction  effects  are  considered  to  be  strong  enough  for 
all  frequencies  so  that  the  reflected  ray  is  sufficiently 
reduced  in  amplitude  to  be  negligible. 

Also  neglected  are  reflections  from  planes  that  are  either  very 
short  or  very  remote  so  that  the  contribution  to  the  sound 
intensity  at  the  receiver  is  less  than  10%  of  the  intensity 
received  via  direct  (or  diffracted)  rays  from  the  road  segment 
under  consideration.  The  analytical  formulation  for  this 
criterion  is 

d’Aa  J 

distance  from  the  receiver  to  the  road  segment, 

aspect  angle  of  the  road  segment  at  the  receiver, 

distance  from  the  road  segment  to  a receiver 
location  imaged  about  the  reflector. 

aspect  angle  of  the  barrier  at  the  image  receiver, 

attenuation  in  dB  by  diffraction  of  the  direct  ray 
due  to  a possible  barrier  (referred  to  a frequency 
of  500  Hz) . 


where  d = 
Aa  = 
d'  = 

Aa ' = 
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Single  reflections  are  considered  exclusively;  contributions  from 
rays  that  strike  two  or  more  reflectors  are  ignored.  It  is  essen- 
tial for  the  calculation  of  higher  order  statistical  parameters 
of  road  traffic  noise  that  the  reflections  of  sound  from  a certain 
road  segment  be  treated  as  amplifications  .of  the  direct  (or  dif- 
fracted) rays  and  not  as  uncorrelated  contributions  from  inde- 
pendent road  segments.  The  factor  F multiplying  the  intensity  of 
the  direct  sound  from  a road  segment  is  calculated  using 


where  the  subscript  i indicates  reflections  at  N different 
surfaces,  each  of  which  might  be  diffracted  by  a barrier  before 
or  after  reflection  and  therefore  might  have  an  attenuation  . 
The  factor  F is  calculated  as  a function  of  frequency.  The 
notation  in  Eq.  (54)  is  the  same  as  in  Eq.  (53)  except  for 
the  angle  Aa*  which  denotes  the  aspect  angle  of  the  road  seg- 
ment at  the  image  receiver  (see  Fig.  9)* 

2.9  Combination  of  Attenuation  and  Reflection 
2.9.1  Atmospheric  absorption 

The  computer  program  accounts  for  atmospheric  absorption  in 
combination  with  barrier  diffraction,  ground  absorption  and 
reflections.  The  path  length  used  for  calculating  the  atmo- 
spheric absorption  is  the  direct  distance  from  the  source  to 
the  receiver  and  is  not  corrected  for  the  path  length  difference 
6 of  diffracted  rays  nor  for  the  increased  path  length  of  re- 
flected rays.  The  factor  F,  defined  by  Eq . (54),  is  multiplied 


F = 10  dB/10  + l 


Aa!  d 

l 


(54) 


i=l  A ad! 

l 
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ROAD  SEGMENT 


RECEIVER 


FIG. 9 EXAMPLE  FOR  RAY  TRACES  FROM  A REFLECTOR, 
INDICATING  THE  CONSTRUCTION  BY  MEANS  OF 
AN  IMAGE  RECEIVER 
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by  10  DA/10,  where  is  defined  by  Eq.  (43), 
culate  for  each  individual  road  segment  s the 

10"D/1°  = F 10-DA/10  10-DG/10  , 

which  is  employed  for  the  calculation  of  the  energy  mean  level 
in  Eqs . (16)  and  (17),  and  of  the  cumulants  in  Eq.  (18). 

2.9.2  Ground  absorption 

The  computer  program  accounts  for  ground  absorption  only  if  the 
rays  travelling  over  the  ground  are  not  diffracted  by  a barrier. 
This  concept  is  based  on  field  experience  with  low  ground  cover 
\_13i14~\.  For  high  ground  cover,  it  might  be  necessary  to  recon- 
sider the  approach  taken.  As  shown  with  the  formula  Eq . (55), 
the  same  ground  absorption  is  assumed  for  the  direct  rays  and 
for  all  reflections. 

2.9.3  Re  flection 

The  computer  program  accounts  for  reflections  in  combination 
with  diffraction  provided  that  there  is  only  a single  diffrac- 
tion before  or  after  the  reflection  and  that  the  path  length 
increase  due  to  diffraction  is  less  than  12.5  ft.  Doubly 
diffracted  reflections  are  neglected  as  well  as  very  weak 
single  reflections  that  suffer,  in  the  500  Hz  band,  the 
maximum  attenuation  of  24  dB  assumed  for  barrier  diffraction. 

The  attenuation  of  reflected  rays  by  diffraction  is  calculated 
for  one  location  on  the  road  segment  only:  the  point  nearest 

to  the  image  receiver.  No  attempt  has  been  made  to  refine 
this  calculation  by  checking  for  the  attenuation  from  other 
points  on  the  road  segment,  since  the  contribution  of  dif- 
fracted reflections  will  be  generally  small  and,  hence, 

>s. 

inaccuracies  of  the  calculation  will  be  negligible. 


in  order  to  cal- 
factor 

(55) 
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3. 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


3.1  Introduction 

The  computer  program  is  designed  to  run  on  a general  purpose 
computer  of  the  IBM  7094  type  such  as  the  one  currently  installed 
in  the  Computer  Technology  Division  at  TSC.  The  program  is 
written  in  FORTRAN-IV  and  is  designed  to  run  in  batch  mode.  Be- 
cause of  the  modular  organization  of  the  program,  only  a modest 
level  of  effort  should  be  required  to  modify  it  for  an  inter- 
active mode  of  operation.  As  currently  implemented,  inputs  are 
provided  via  punched  cards  and  outputs  are  provided  on  a line 
printer  or  any  other  device  specified. 

The  organizational  scheme  of  the  program  is  shown  in  the  follow- 
ing flow  diagram.  The  major  blocks  of  code  are  contained  in  the 
main  program  and  in  the  subroutines  INPUT  and  GEOMRY . The  sub- 
routine INPUT  reads  input  data  provided  on  cards,  stores  this 
information  in  core  memory,  and  outputs  it  on  a line  printer. 

The  subroutine  GEOMRY  provides  the  bulk  of  the  calculations  per- 
formed by  the  computer  program  and  uses  a number  of  additional 
subroutines  and  functions  in  the  process.  Additional  subroutines 
and  functions  called  by  the  main  program  are  FILES,  which  assigns 
device  numbers  to  the  various  I/O  devices,  INTER,  which  computes 
reference  intensities  for  t'he  various  vehicle  groups  specified  by 
the  user,  XPROB,  which  serves  for  a refined  calculation  of  the 
percentile  level  L^q,  and  GAUSS,  which  checks  the  calculated 
percentile  levels  against  the  background  noise. 

The  main  program  and  the  subroutines  INPUT  and  GEOMRY  are  dis- 
cussed individually  in  some  detail  in  the  remainder  of  this  section 
and  the  remaining  subroutines  and  functions  are  described  very 
briefly.  Listings  for  the  entire  computer  program  are  given  in 
Sec.  3.6.  Details  on  program  usage  are  provided  in  Sec.  4. 
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ORGANIZATIONAL  SCHEME  OF  COMPUTER  PROGRAM 


FILES 


Main  Program 


INPUT 

INTER 

GEOMRY 

GAUSS 

Subrout ines : 

Functions : 

BLOCN 

ACOS 

ENDPT 

AM  AG 

IMAGE 

ANGLE 

INTCPT 

BARFAC 

MOVE 

DEL 

NRPT 

DSQR 

NRI 

HEIGHT 

REPLCE 

IEPS 

SECTN 

KCUT 

TRI 

KPOS 

TAN 

ZCOR 

XPROB 

XINGAM 

GAMMA 


! 


The  following  definitions  are  used  throughout  the  remainder  of 
this  manual.  A "road"  or  "roadway"  is  defined  as  a set  of 
connected  straight  "road  sections"  for  which  the  traffic  condi- 
tions are  uniform.  That  is,  the  number  of  vehicle  groups  and 
the  concentration  and  mean  vehicle  speed  associated  with  each 
group  are  the  same  for  all  sections  of  a given  road.  Barriers 
are  treated  similarly.  Thus,  a "barrier"  is  defined  as  a set 
of  connected  straight  "barrier  sections"  for  which  the  absorptive/ 
reflective  properties  are  the  same. 

3.2  Main  Program 
3.2.1  Overview 

The  primary  function  of  the  main  program  is  to  control  the  flow 
of  operations.  Various  subroutines  are  called  upon  to  perform 
the  bulk  of  the  input-output  operations  and  the  analysis.  In 
addition,  certain  variables  are  initialized  in  the  main  program 
and  the  final  computation  of  energy  mean  level,  of  the  noise 
pollution  level,  and  of  percentile  levels  of  the  A-weighted  SPL 
is  made. 

The  final  computation  starts  with  the  combination  of  cumulants 
of  the  sound  intensity  contributed  from  the  traffic  under  consid- 
eration and  from  the  background  noise  input  for  each  receiver. 

For  the  subsequent  calculation  of  statistical  parameters,  the 
higher  order  cumulants  are  normalized  with  the  first  cumulant, 
which  represents  the  total  energy  mean. 

Following  the  calculation  of  the  noise  pollution  level  and  of 
percentile  levels  with  Eqs.  (4)  through  (8)  of  Sec.  2.1,  a refined 
calculation  of  L is  performed  with  the  evaluation  of  the  cumu- 
lative distribution  given  in  Eq.  (39)  of  Sec.  2.4.  Probabilities 
are  calculated  as  a function  of  the  sound  pressure  level,  starting 
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Main  Program 
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Main  Program  (Continued) 


41 


with  the  energy  mean  level  Lg  and  descending  in  steps  of  1 dB  if 
the  probability  for  Lg  exceeds  0.9,  or  ascending  in  steps  of  1 dB 
if  the  probability  for  Lg  is  less  than  0.9.  The  result  for  the 
level  L10  follows  from  a cross-over  of  the  step-wise  calculated 
probability  and  the  probability  0.9. 

When  a background  noise  level  is  specified,  the  main  program 
finally  calculates  the  increase  of  all  calculated  levels  above 
the  normally  distributed  levels  of  background  noise. 

A flow  diagram  of  the  main  program  is  given  below.  This  diagram 
is  self-explanatory.  Whenever  possible,  statement  numbers  are 
shown  so  that  reference  may  be  made  to  the  corresponding  block 
of  code  in  the  program  listing. 

3.2.2  List  of  variables  for  main  program 

Variables  appearing  in  the  main  program  are  listed  below.  Array 
variables  are  not  generally  identified  as  such;  the  dimension- 
ality of  the  variables  can  be  determined  from  the  program 
listing  in  Sec.  3-6.  The  variables  XR1  and  XR2  listed  below 
are  vector  quantities  which  indicate  the  coordinates  of  road 
end  points.  The  first  element  of  the  vector  indicates  the 
x-coordinate , the  second  element  the  y-coordinate , and  the 
third  element  the  z-coordinate . This  format  is  used  widely 
throughout  the  remainder  of  the  entire  computer  program  to 
denote  coordinates. 

A1  Temporary  variable 

A2  Temporary  variable 

AA3  Temporary  variable 

AA4 


Temporary  variable 


accounting  for  standard  deviation  of  reference 


AI  Index 

CO  Factor 
level 

FLAG  Indicator  for  L above  or  below  energy  mean  level 
I Index 


IBAR  Barrier  number 


DXL10 

Increase 

in 

level 

Li» 

above 

background 

noise 

DXL50 

Increase 

in 

level 

Lo 

above 

background 

noise 

DXL90 

Increase 

in 

level 

L,0 

above 

background 

noise 

DXLA 

Increase 

in 

energy 

mean  level  above  background  noise 

DXLNP 

Increase 

in 

noise 

pollution 

level  above 

■ background  noise 

IERR  Error  index 
II  Index 
IQ  Index 


ISEG  Barrier  section  number 
J Index 

M Roadway  number 
N Road  section  number 


NB 

Number 

of 

barriers 

NF 

Number 

of 

frequency  bands 

NG 

Number 

of 

absorptive  ground 

strips 

NLIM 

Number 

of 

sections  defining 

a roadway 

NQ 

Number 

of 

vehicle  types 

NQQ 

Number 

of 

groups  within  one 

vehicle  type 
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NQS 

NR 

NRC 

NRSM1 

PRO 

PR 

R 

RDIN 

RX 

RY 

RZ 

S 

SIG 

SIGL 

STB 

TEMPI 

TEMP2 

VEXPH 

XI 

XKA 

XKAP 

XLA 

XLB 


Vector  notation  for  number  of  vehicle  groups 

Number  of  roadways 

Number  of  receivers 

Number  of  sections  for  one  roadway 

Probability  of  total  noise 

Probability  of  background  noise 

Parameter  defined  in  Eq.  (36) 

Vector  notation  for  initialization  parameters 
x-coordinate  of  roadway  point 
y-coordinate  of  roadway  point 
z-coordinate  of  roadway  point 
Parameter  defined  in  Eq.  (37) 

Standard  deviation  of  reference  level 

Standard  deviation  of  A-weighted  sound  pressure  level 
Standard  deviation  of  A-weighted  background  noise  level 
Temporary  variable 
Temporary  variable 
Vehicles  per  foot 

Sound  Intensity  normalized  by  its  mean 
Cumulant  of  the  A-weighted  sound  intensity 
Normalized  cumulant  of  the  A-weighted  sound  intensity 
A-weighted  overall  level 

Mean  value  of  A-weighted  background  noise  level 
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XLE 

A-weighted 

intensity  and 

level 

in 

frequency  bands 

XL10 

A-weighted 

SPL 

exceeded 

10% 

of 

the 

time 

XL50 

A-weighted 

SPL 

exceeded 

50% 

of 

the 

time 

XL90 

A-weighted 

SPL 

exceeded 

90% 

of 

the 

time 

XLNP  Noise  pollution  level 

XMPH  Vehicle  speed  in  mph 

XR1  Road  section  initial  point 

XR2  Road  section  end  point 

XRC  x-coordinate  of  receiver 

YRC  y-coordinate  of  receiver 

Z1  Normalized  sound  intensity 

ZI  A-weighted  sound  pressure  level 

ZRC  z-coordinate  of  receiver 

3. 3 Subrouti ne  INPUT 
3.3.1  Overview 

This  subroutine  reads  the  input  data  provided  by  the  user  on 
punched  cards,  stores  these  data  in  core,  and  outputs  the  data 
on  the  user  teletype  so  that  a permanent  record  is  created. 

The  computer  program  is  organized  in  such  a way  that  the  user 
need  not  redefine  the  entire  situation  each  time  a new  case  is 
analyzed.  Instead,  it  is  necessary  to  specify  only  those 
aspects  of  the  problem  that  are  different  from  the  previous 
case.  Accordingly,  each  block  of  input  data  is  preceded  by  a 
control  card  that  indicates  the  type  of  data  (e.g.,  barrier 


specifications)  and  the  number  of  items  (e.g.,  number  of  barriers) 
to  follow.  After  a special  control  character  has  signified  the 
end  of  data  input,  the  entire  set  of  input  variables  is  typed 
out  and  control  returns  to  the  main  program. 

A flow  diagram  of  this  subroutine  is  given  below,  followed  by  a 
list  of  variables. 


HS 


Subroutine  Input 
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For  IGO 


3: 


For 


For 


IGO 


IGO 
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For  IGO  = 6 


50  <> 

For  all  roadways: 
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3.3.2  List  of  variables  for  subroutine  INPUT 


ALPHA 

Alphanumeric  information 

BGS 

Width  of  absorptive  ground  strip 

BX 

x-coordinate  of  barrier  point 

BY 

y-coordinate  of  barrier  point 

BZ 

z-coordinate  of  barrier  point 

I 

Index 

/ 

II 

Number  of  items  in  data  block 

12 

Dummy  variable 

IA 

Alphanumeric  "A" 

IB LAST 

Barrier  type 

IDN 

Index  for  program  initialization  parameter 

IDUM 

Index  for  type  of  absorptive  ground  cover 

IG 

Alphanumeric  "G" 

IGO 

Index  for  data  blocks 

ILAST 

Alphanumeric  indicator  for  last  section 

IR 

Alphanumeric  "R" 

IT 

Alphanumeric  nT" 

ITY 

Vehicle  type 

J 

Index 

K 

Index 

LAST 

Alphanumeric  "L" 

NB 

Number  of  barriers 
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NBSM1 

NF 

NG 

NQ 

NQC 

NQS 

NR 

NRC 

NRSM1 

NSEC 

NSM1 

NTEMP 

RDIN 

RX 

RY 

RZ 

STB 

VALUE 

VEH 

VEXPH 

VTEMP 

XLB 


j 

i 

! 

! 

I 


Number  of  sections  for  one  barrier 
Number  of  frequency  bands 
Number  of  absorptive  ground  strips 
Number  of  vehicle  types 

Number  of  groups  within  one  vehicle  type 

Vector  notation  for  number  of  vehicle  groups 

Number  of  roadways 

Number  of  receivers 

Number  of  sections  for  one  roadway 

Section  number 

NSEC-1 

Temporary  variable 

Vector  notation  for  initialization  parameters 
x-coordinate  of  roadway  point 
y-coordinate  of  roadway  point 
z-coordinate  of  roadway  point 

Standard  deviation  of  A-weighted  background  noise  level 
Initialization  parameter 
Vehicles  per  hour 

Vehicles  per  foot,  vehicles  per  hour 
Vehicle  speed 

Mean  value  of  A-weighted  background  noise  level 
Speed  in  mph  for  one  group  of  vehicles 
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XMPH 

Vector  notation  for  vehicle  speed 

XNIGHT 

Height  adjus 

tment  for  receiver  coordinate 

XRC 

x-coordinate 

of  receiver 

XXG1 

x-coordinate 

for  point  of 

ground 

strip 

center 

line 

YRC 

y-coordinat e 

of  receiver 

YYG1 

y-coordinate 

for  point  of 

ground 

strip 

center 

line 

ZRC 

z-coordinate 

of  receiver 

ZZG1 

z-c  oordinat e 

for  point  of 

ground 

strip 

center 

line 

3.4  Subroutine  GEOMRY 
3.4.1  Overview 

The  bulk  of  the  computational  effort  of  this  subroutine  Is 
devoted  to  finding  segments  of  a straight  road  section  for 
which  sound  from  all  points  along  the  segment  reaching  a par- 
ticular receiver  is  nearly  equally  attenuated  by  barrier  dif- 
fraction or  ground  absorption  or  is  equally  amplified  by 
reflections . 

Because  of  the  complexity  of  this  subroutine,  a detailed  de- 
scription with  a series  of  flow  diagrams  is  presented  in 
Sec.  3.^.2.  Before  reviewing  these  diagrams,  however,  the 
reader  should  become  familiarized  with  the  basic  philosophy 
of  the  computational  scheme  outlined  below. 

As  shown  in  the  following  diagram  giving  an  overview  of  sub- 
routine GEOMRY,  checks  are  performed  on  two  different  levels  — 
first  on  all  the  barriers,  then  on  all  the  absorptive  ground 
strips,  and  finally  on  all  the  reflectors.  On  the  lower  level, 
preliminary  checks  are  made  mostly  in  the  two  dimensions  of 
the  x-y  plane  to  determine  whether  a barrier,  a ground  strip. 
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AN  OVERVIEW  OF  SUBROUTINE  GEOMRY 


I 
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or  a reflector  can  be  at  all  effective  for  the  road  segment 
and  the  receiver  under  consideration.  During  these  checks, 
arrays  are  created  for  a smaller  number  of  barriers,  ground 
strips,  and  reflectors  that  need  to  be  investigated  more 
closely  at  the  higher  level. 

At  this  higher  level,  break  points  between  road  segments  pos- 
sibly affected  and  unaffected  by  barriers,  ground  strips,  or 
reflectors  are  determined  from  checks  first  made  in  the  two 
dimensions  of  the  x-y  plane.  If  the  road  segment  between 
the  initial  point  and  the  break  point  is  unaffected,  then 
no  further  checks  in  three  dimensions  are  made  on  the  actual 
effect  of  the  remaining  road  segment,  but  the  break  point  is 
considered  as  the  new  end  point  of  the  road  segment.  The 
remaining  road  segment  is  investigated  after  the  calculations 
for  the  first  road  segment  are  finished. 

If  the  road  segment  between  the  initial  point  and  the  break 
point  is  possibly  affected  by  a barrier,  a ground  strip,  or 
a reflector,  the  break  point  is  taken  as  a preliminary  end 
point  until  the  subsequent  analysis  in  three  dimension  shows 
whether  there  is  an  actual  effect  that  makes  it  necessary  to 
keep  the  break  point  as  an  end  point  or  whether  there  is  no 
effect  and  the  break  point  can  be  dropped. 

Since  all  the  checks  are  performed  first  for  diffraction,  then 
for  absorptive  ground  strips,  and  finally  for  reflections,  the 
initial  road  section  may  be  reduced  several  times  by  each  group 
of  elements.  The  procedure  of  subdividing  a road  section  is 
diagrammed  in  Pig.  10.  The  final  reduction  in  each  group  is 
stored  in  terms  of  end  points  XR2D  for  barriers,  XR2G  for 
absorptive  ground  strips,  and  XR2  for  reflectors.  By  this 
means,  duplicate  calculations  are  avoided.  After  the  noise 
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CONSTANT  ATTENUATION  BY  DIFFRACTION 
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CONSTANT  ATTENUATION  BY  GROUND  ABSORPTION 
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UNIFORM  CONTRIBUTION  OF  REFLECTED  RAYS 
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XR2 


XR2G 


XR2D 


XR20 


FIG. 10  SUBDIVISION  OF  ROAD  SECTION 


from  the  road  segment  XR10,  XR2  is  analyzed,  the  computation 
continues  for  possible  reflections  from  the  road  segment  XR2 , 
XR2G.  After  the  noise  from  the  road  segment  XR1 , XR2G  is 
completely  analyzed,  the  computation  continues  with  considera- 

i 

tion  of  ground  absorption  and  reflections  for  the  road  segment 
XR2G,  XR2D.  The  designation  of  end  points  during  the  course 
of  this  analysis  is  illustrated  in  Pig.  11.  Finally,  the 
computation  continues  with  the  entire  problem  of  barrier  dif- 
fraction, ground  absorption,  and  reflection  for  the  remaining 
road  segment  XR2D,  XR20  until  the  noise  from  the  entire  road 
section  XR10,  XR20  has  been  analyzed.  At  that  point,  program 
control  is  returned  to  the  main  program. 

3.4.2  Detailed  flow  diagrams 

The  flow  diagrams  presented  at  the  end  of  this  section  show  the 
operations  and  important  branch  points  of  the  computations  that 
are  indicated  in  the  overview  diagram  presented  above.  Again, 
statement  numbers  are  shown  so  that  reference  may  be  made  to 
the  corresponding  block  of  code  in  the  program  listing  given 
in  Sec.  3*6.  The  following  comments  are  provided  to  supplement 
the  information  contained  in  the  flow  diagrams. 

Preliminary  computations 

The  main  computation  shown  in  the  first  flow  chart  concerns  the 
distance  between  the  near  point  XNPT  on  the  source  line  and  the 
receiver  XRC.  This  distance  enters  into  the  formula  for  the 
mean  energy  level  and  is  a quantity  that  is  characteristic  of 
the  entire  road  section. 
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CONSTANT  ATTENUATION  BY  GROUND  ABSORPTION 


UNIFORM  CONTRIBUTION  OF  REFLECTED  RAYS 
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FIG.  11 


FURTHER  SUBDIVISION  OF  ROAD  SECTION 


Preliminary  selection  of  diffracting  barriers 


The  computations  are  made  in  the  two  dimensions  of  the  x-y 
plane,  and  no  checks  are  made  for  the  height  of  barriers.  A 
barrier  whose  top  line  crosses  the  roadway  in  a projection  on 
the  x-y  plane  is  considered  illegal  and  causes  the  program  to 
abort.  Thus,  even  a barrier  crossing  underneath  an  elevated 
roadway  is  illegal.  A barrier  stored  in  the  array  B1  is  not 
necessarily  high  enough  actually  to  diffract  sound  rays 
travelling  from  the  road  section  to  the  receiver.  The  storage 
of  the  relative  location  (see  SUBROUTINE  BLOCN)  saves  consider- 
able computation  for  very  high  barriers  that  shield  the  entire 
road  section  against  the  receiver. 


Preliminary  selection  of  ab sorbing  ground  strips 

The  computations  are  entirely  analogous  to  the  preliminary 
selection  of  diffracting  barriers.  No  check  is  made  for  the 
distance  of  the  ground  strip  center  line  from  the  road  section. 
Therefore,  a road  section,  or  a part  of  it,  which  is  inadver- 
tently specified  on  a ground  strip  will  not  be  detected  in  all 
cases  by  the  program. 

Diffraction  of  direct  ray 

The  initialization  of  a path  length  DELP0  = -0.2  means  that  the 
attenuation  by  diffraction  is  set  equal  to  zero  in  the  frequency 
band  centered  at  500  Hz. 

The  test  for  whether  or  not  a barrier  possibly  shields  the 
entire  road  is  made  in  the  x-y  plane  and  does  not  account  for 
the  actual  height  of  the  barrier. 
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barrier 


The  height  difference  for  sound  rays  from  cars  above  the 
provides  a rather  conservative  check  for  possible  effects  of 
diffraction.  A positive  height  difference  means  the  source 
is  visible  from  the  receiver  location. 

Only  the  most  effective  of  a number  of  parallel  barriers  is 
considered;  thus,  the  barrier  having  the  strongest  diffraction 
is  retained.  Note  that  this  barrier  is  not  necessarily  the 
highest  one.  It  is  just  as  likely  to  be  the  barrier  nearest 
to  the  road  or  to  the  receiver. 


The  calculation  of  a road  segment  with  sufficiently  similar 
diffraction  is  performed  for  cars  only.  It  is  then  assumed 
that  the  same  result  holds  for  trucks  having  a different 
source  height. 


The  check  for  maximum  diffraction  of  sound  in  the  frequency 
band  centered  at  500  Hz  is  made  for  trucks  only,  because  the 
higher  source  height  of  trucks  will  result  in  smaller  dif- 
fraction than  for  cars.  If  this  test  reveals  a path  length 
difference  of  more  than  12.5  ft,  further  checks  on  other 
barriers  can  be  skipped  for  the  road  segment  considered. 

If  a very  high  barrier  shields  the  entire  remaining  road 
section,  no  further  tests  for  diffraction  of  the  direct  ray 
are  made.  This  procedure  may  result  in  some  inaccuracies 
for  very  long  road  sections  for  which  the  diffraction  of 
sound  from  points  near  to  the  receiver  is  very  strong  but 
for  which  the  diffraction  of  sound  from  remote  points  might 
be  considerably  weaker.  However,  these  inaccuracies  are  no 
more  serious  than  those  resulting  from  the  assumption  of  the 
arbitrary  value  of  12.5  ft  for  the  maximum  path  length  dif- 
ference, and  are  thus  neglected. 
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Ground,  ab sorption 

Absorptive  ground  strips  are  disregarded  for  sound  from  a road 
segment  whenever  an  attenuation  different  from  zero  is  found 
due  to  barrier  diffraction.  The  criterion  for  this  decision 
may  be  modified  in  the  future  to  be  5 or  10  dB , corresponding 
to  path  length  differences  of  0 or  0.1  ft  at  500  Hz,  if  field 
experience  shows  that  this  is  necessary.  Furthermore,  the 
decision  to  neglect  ground  strips,  which  is  related  to  certain 
heights  of  rays  above  the  center  line  of  the  ground  strip, 
is  not  based  on  experimental  evidence  and  may  require  modi- 
fications . 

Preliminary  selection  of  reflectors 

By  the  introduction  of  an  image  receiver  location,  the  reflec- 
tion problem  becomes  similar  to  the  diffraction  problem.  A 
reflector  in  the  path  from  the  road  segment  to  an  image  receiver 
is  effective  whenever  a barrier  in  the  path  from  the  road  seg- 
ment to  a receiver  strongly  diffracts  the  sound. 

In  addition  to  the  preliminary  checks  made  above  for  diffracting 
barriers,  the  intensity  of  the  direct  (or  diffracted)  sound 
from  the  road  segment  considered  is  compared  to  the  potential 
maximum  contribution  of  each  reflection.  This  check  assures 
that  only  essential  reflectors  are  considered.  Since  a re- 
flector might  not  be  high  enough  or  reflections  might  be 
strongly  attenuated  by  diffraction  at  additional  barriers 
in  the  ray  path,  the  reflectors  found  at  this  stage  are  con- 
sidered "potential"  reflectors. 
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Also  involved  in  the  part  entitled  "Preliminary  Selection  of 
Reflectors"  is  the  preliminary  selection  of  barriers  which 
can  possibly  diffract  the  reflected  sound.  The  check  per- 
formed is  based  on  the  length  and  orientation  of  a barrier 
section  and  on  the  distance  from  the  image  receiver  relative 
to  the  respective  parameters  of  the  road  section  and  the 
receiver  considered. 

Calculation  of  potential  reflections 

The  computations  in  this  section  are  analogous  to  those  made 
for  the  diffraction  problem,  except  for  a simplified  decision 
concerning  the  height  of  the  reflector.  If  the  height  dif- 
ference above  the  reflector  of  a sound  ray  from  the  nearest 
point  on  the  road  segment  leads  to  the  acceptance  or  rejection 
of  the  reflector,  no  further  check  is  made  for  other  source 
points  on  the  road  segment. 

Since  diffraction  of  the  reflected  ray  is  not  yet  checked,  a 
reflection  is  still  called  a potential  one  and  the  end  point 
of  the  road  segment  is  preliminary. 

Consideration  of  diffraction  before  reflection 

Checks  for  barriers  in  the  area  defined  by  projections  of  the 
four  points  XR1,  XR2 , XRB3 , and  XRB4  into  the  x-y  plane  are 
made  by  considering  first  the  triangle  containing  image  rays 
and  then  the  triangle  formed  by  the  road  segment  and  the  image 
receiver . 
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If  a barrier  is  found  to  be  high  enough  for  possible  diffraction, 
the  reflector  is  checked  again  to  determine  whether  or  not  it 
is  high  enough  to  reflect  the  diffracted  rays  (which  now  come 
from  an  effective  source  that  might  be  considerably  higher 
than  the  roadway). 

Calculation  of  the  path  length  difference  between  diffracted 
and  direct  rays  from  only  the  near  point  of  the  road  segment 
implies  that  the  diffraction  of  sound  rays  from  other  source 
points  is  about  the  same. 

Very  strongly  diffracted  reflections  are  neglected.  The  deci- 
sion is  made  on  the  basis  of  the  diffraction  of  sound  from 
trucks,  since  rays  from  cars  are  even  more  strongly  diffracted. 

Consideration  of  diffraction  after  reflection 

Checks  are  made  for  diffracting  barriers  in  the  triangle  defined 
by  projections  of  the  reflector  segment  XRB3,  XRB4  and  the  re- 
ceiver XRC  onto  the  x-y  plane. 

After  a barrier  has  been  found  which  is  high  enough  for  possible 
diffraction,  the  reflector  is  checked  to  determine  whether  or 
not  it  is  high  enough  to  reflect  sound  towards  the  top  line 
of  the  diffracting  barrier,  which  might  be  considerably  higher 
than  the  receiver. 

The  calculation  of  the  path  length  difference  implies  simplify- 
ing assumptions  similar  to  those  for  the  problem  of  diffraction 
before  reflection. 

Reflections  are  neglected  in  the  case  of  diffraction  before  and 
after  reflection  and  in  the  case  of  very  strong  diffraction 
after  reflection. 
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Final  computations  and  summation  of  oumulants  of  the  sound  intensity 

Since  contributions  from  reflections  to  the  sound  intensity 
received  from  one  road  segment  are  not  statistically  independent, 
they  are  not  entered  directly  in  the  calculation  of  cumulants, 
but  they  are  used  to  calculate  a gain  factor  of  the  direct  sound 
from  a road  segment.  The  subsequent  computations  are  directly 
oriented  at  the  formulas  of  Eq.  (1 6)  for  the  mean  sound  intensity 
and  of  Eq.  (18)  for  the  higher  order  cumulants  of  the  sound 
intensity . 

After  the  calculation  for  a road  segment  has  been  completed,  a 
point  1 ft  beyond  the  end  point  of  the  segment  is  taken  as  the 
initial  point  of  the  next  road  segment.  This  procedure  avoids 
computational  complications  without  unduly  sacrificing  accuracy. 
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Preliminary  Computations 
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Preliminary  Selection  of  Diffracting  Barriers 


For  all  sections  of  all  barriers: 


65 


Preliminary  Selection  of  Absorbing  Ground  Strips 


For  all  ground  strips : 


Does  center  line  of  the  ground  \ Y 
strip  cross  the  road  section? 


-Q  RETURN  ^ 


Does  the  center  line  of  the 
ground  strip  or  a part  of  it 
lie  in  the  triangle  formed  by 
the  road  section  XR10,  XR20 
and  the  receiver  XRC? 


Store  the  information  about 

the  ground  strip  and  the 
relative  location  in  the 
arrays  TAI . BGT,  and  KGCODE 

Disregard  the 
ground 
strip 

1 8 
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Diffraction  of  Direct  Ray 


For  all  vehicle  types: 


Initialize  path  length  difference  between  a dif- 
fracted ray  and  a direct  ray  to  DELP0  = -0.2 


For  all  barriers  in 
array  Bl: 


30  lh 


4i_j < Does  barrier  possibly  shield  the 

* entire  road  section  XR10,  XR20? 


Find  break  point  XK  between  possibly 
shielded  and  unshielded  road  segments 


Use  road  segment 
endpoint  XR2  as 
preliminary  end- 
point XK 


Reduce  road  seg- 
ment considered 
to  XR1 , XR2  = XK 


Is  initial  segment  XR1 , 
XK  possibly  shielded? 
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,Y 


For  all  vehicle  types : 


Find  height  HDIFF  above  barrier  of  sound 
ray  from  nearest  point  on  the  road 
segment  XN1  to  receiver  XRC 


For  cars  only : 
Is  HDIFF  > 20? 


Disregard  barr. 

ier 
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N 


50 


Find  path  length  difference  DELP 
between  diffracted  and  direct 
ray  from  XN1  and  XRC 


N 


Is  diffraction  stronger  than 
previously  calculated? 


Has  barrier  been  found  most 
important  for  cars? 
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67 


68 


Ground  Absorption 
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Preliminary  Selection  of  Reflectors 

For  all  sections  of  all  barriers: 
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Calculation  of  Potential  Reflections 
For  all  reflectors  in  the  array  Rl: 
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Consideration  of  Diffraction  Before  Reflection 
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Consideration  of  Diffraction  After  Reflection 
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Final  Computations  and  Summation  of  Cumulants  of  the  Sound  Intensity 
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24  6 

Initialize  for 

remaining  road 

segments:  ANG1  = |ANG2| 

i 

Is  endpoint  XR2  of  road  seg- 
ment within  1 ft  of  endpoint 
XR20  of  road  section? 


N 


Move  endpoint  XR1  towards 
XR20  1 ft  beyond  XR2 


Is  endpoint  XR2  within 
1 ft  of  XR2G? 


9 1 


2 4 7 0 


N 


Continue  with 
reflection 
problem  for 
remaining 
road  segment 


Continue  with 
ground  absorp- 

Is  endpoint  XR2  within  >\v72 

1 ft  of  XR2D? 

tion  problem 
for  remaining 
road  segment 

25< 

,N 

Continue  with  diffraction  of  direct 
ray  for  remaining  road  segment 
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3.4.3  List  of  variables  for  subroutine  GEOMRY 


A 

Atmospheric  attenuation  factor,  ground  absorption 
parameter 

ADST 

ANG/DIST 

ANG 

Angle  subtended  at  receiver  by  road  segment 

ANG1 

Angle  a in  Pig.  1 

ANG2 

Angle  a2  in  Pig.  1 

ANGI 

Angle  subtended  at  image  receiver  by  road  segment 

ANGIMG 

Angle  subtended  at  image  receiver  by  barrier  section 

B1 

Array  of  potentially  diffracting  barrier  segments 

BGS 

Width  of  absorptive  ground  strip 

BGT 

Width  of  selected  absorptive  ground  strip 

BX 

x-coordinate  of  barrier  point 

BY 

y-coordinate  of  barrier  point 

BZ 

z-coordinate  of  barrier  point 

CA1 

cos  ( a x ) 

CA2 

cos  (a  ) 
2 

CPREV 

Angular  function  C defined  in  Eq.  (19) 

CQ 

Factor  accounting  for  standard  deviation  of  reference 

DELM 

Mean  path  length  difference 

DELP 

Path  length  difference 

DELP0 

Maximum  path  length  difference  for  diffraction  of 
direct  ray 

DELP1 

Maximum  path  length  difference  for  diffraction  before 
reflection 
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DELP2 

DELPA 

DELR 

DELTA 

DIST 

DISTI 

DISTJ 

DL 

DN1 

DN1I 

DN2 

DR1 

DRK 

PB 

PCTR 

PG 

HDIPA 


Maximum  path  length  difference  for  diffraction  after 
reflection 

Path  length  difference  for  diffraction  of  direct  ray 
Path  length  difference  for  reflected  ray 
Distance  along  the  roadway 

Distance  from  the  receiver  to  the  source  line 

Distance  from  the  image  receiver  to  the  source  line 

Distance  from  the  image  receiver  to  the  diffracting 
barrier 

Mean  path  length  over  an  absorptive  ground  strip 

Distance  from  the  receiver  to  the  nearest  point  of 
the  road  segment 

Distance  from  the  image  receiver  to  the  nearest  point 
of  the  road  segment 

Distance  from  the  image  receiver  to  the  nearest  point 
of  the  diffracting  barrier 

Distance  from  the  receiver  to  the  initial  point  of 
the  road  segment 

Distance  from  the  receiver  to  the  preliminary  end  point 
of  the  road  segment 

Attenuation  factor  accounting  for  diffraction  and 
reflections 

Weighting  factor  for  reflections 
Ground  attenuation  factor 

Height  of  ray  from  source  point  XR1  to  receiver  XRC 
above  barrier 
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HDIFF 

Height  of  ray  above  barrier,  reflector,  or  ground 
strip 

HGA 

Data  for  effective  height  of  ground  cover 

I 

Index 

IA 

Alphanumeric  "A" 

IBAR 

Barrier  number 

IBLAST 

Barrier  type 

ICODE 

Number  for  intermediate  print  out 

IDUM 

Index  for  kind  of  absorptive  ground  cover 

IDXR 

Number  of  reflections 

IERR 

Error  index 

IGRA 

Ground  strip  number 

IF 

Index  for  frequency  bands 

II 

Index 

IK 

Frequency  band  number 

IKIN 

Index  "for  kind  of  absorptive  ground  cover 

IP 

Frequency  band  number 

IQ 

Index 

IR 

Alphanumeric  "R" 

ISEG 

Barrier  section  number 

IT3 

2N-1 

ITRIG 

Trigger 

KAR 

Alphanumeric  indicator  for  type  of  barrier 
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KBAR 

Total  number  of  barrier  sections;  index  for  barriers 

KBCODE 

Indicator  for  relative  location  of  barrier 

KBAR1 

Reflector  number  in  array  of  selected  reflectors 

KBAR2 

Diffractor  number  in  array  of  barrier  selected  for 

reflected  rays 

KCD 

Indicator  for  relative  location  of  barrier 

KDIFF 

Barrier  number  in  array  of  selected  diffractors 

KP 

Index  for  frequency  bands 

KGA 

Number  of  selected  ground  strips 

KGCODE 

Indicator  for  relative  location  of  ground  strip 

KIMG 

Reflection  number 

KNUMB 

Total  number  of  selected  barrier  sections 

KRDNUM 

Total  number  of  selected  barrier  sections  relevant 

to  reflection 


KREP 

Reflector  number  in  array  of  selected  reflectors 

KRFDF 

Barrier  number  in  array  of  barrier  selected  for 

reflected  rays 

KRNUM 

Total  number  of  selected  reflector  sections 

KTRIG 

Indicator  for  intersection  of  barrier  or  ground  strip 

LOC 

Indicator  for  relative  location  of  barrier  or  ground 

strip 

MDIPF 

Indicator  for  diffraction  before  reflection 

MODD 

Indicator  for  diffraction  of  direct  ray 

MR 

Roadway  number 
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N 


Cumulant  number 


NB 

Number 

of 

NBSEC 

Number 

of 

NBSM1 

Number 

of 

NDIFP 

Number 

of 

NF 

Number 

of 

NG 

Number 

of 

NIMG 

Number 

of 

NLIM 

Number 

of 

NQ 

Number 

of 

NQQ 

Number 

of 

NQS 

Vector 

no' 

NR 

Number 

of 

NREF 

Number 

of 

NRFDF 

Number 

of 

refected  ] 

PP 

Frequency 

R1 

End  point 

RATIO 

Weighting 

RBI 

End  point 

RDIN 

Vector 

no' 

TA1 

End  point 

strip 
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SA1  sin  (a1) 

SA2  sin  (a2) 

T1  Temporary  variable 

T2  Temporary  variable 

T3  Temporary  variable 

VEXPH  Vehicles  per  foot 

VTEMP  Vehicle  speed 

XB1  Initial  point  of  selected  reflector 

XB2  End  point  of  selected  reflector 

XDB1  Initial  point  of  barrier  selected  for  reflected  rays 

XDB2  End  point  of  barrier  selected  for  reflected  rays 

XDB3  Initial  point  of  effective  barrier  segment 

XDB4  End  point  of  effective  barrier  segment 

XG1  Initial  point  of  center  line  of  selected  absorptive 

ground  strip 

XG2  End  point  of  center  line  of  selected  absorptive  ground 
strip 

XG3  Initial  point  of  effective  ground  strip  segment 
XG4  End  point  of  effective  ground  strip  segment 

XIMG  Vector  of  image  receivers  for  all  reflections 

XJ  Preliminary  end  point  of  effective  reflector  segment 

XK  Preliminary  end  point  of  road  segment 

XKA  Cumulant  of  the  A-weighted  sound  intensity 
XKI  Preliminary  end  point  of  image  road  segment 
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XLA 

A-weighted  intensity  in  frequency  bands 

XLE 

Mean  intensity 

XLREF 

Vector  notation  for  reference  intensities 

XN1 

Point  on  road  segment  nearest  to  receiver 

XN1I 

Point  of  road  segment  nearest  to  image  receiver; 
point  on  image  road  segment  nearest  receiver 

XN2 

Point  on  barrier  segment  nearest  to  image  receiver 

XNPT 

Point  on  source  line  nearest  to  receiver 

XNPTI 

Point  on  source  line  nearest  to  image  receiver; 
point  on  image  source  line  nearest  receiver 

XR 

x-coordinate  of  receiver 

XRl 

Initial  point  of  road  segment 

XRLJ2 

Initial  point  of  road  section 

XR1I 

Initial  point  of  image  road  segment 

XR2 

End  point  of  road  segment 

XR2J3 

End  point  of  road  section 

XR2D 

End  point  of  road  segment  with  constant  attenuation 
by  diffraction 

XR2G 

End  point  of  road  segment  with  constant  attenuation 
by  ground  absorption 

XR2I 

End  point  of  image  road  segment 

XRB1 

Initial  point  of  reflector  on  subfile  IREF 

XRB2 

End  point  of  reflector  on  subfile  IREF 

XRB3 

Initial  point  of  effective  reflector  segment 

XRB4 

End  point  of  effective  reflector  segment 

! 
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XRC 


Receiver  point 
XRCI  Image  receiver  point 

XXG1  x-coordinate  of  point  on  ground  strip  center  line 
YR  y-coordinate  of  receiver 
YYG1  y-coordinate  of  point  on  ground  strip  center  line 
ZN10  z-coordinate  of  XN1  or  XN1I 

ZR  z-coordinate  of  receiver 

ZS  Height  adjustment  for  vehicles 
ZZG1  z-coordinate  of  point  on  ground  strip  center  line 


3.5  Description  of  Other  Subroutines  and  Functions 
SUBROUTINE  FILES 

FILES  defines  file  device  numbers  for  the  I/O  devices. 
SUBROUTINE  INTER  (IQ) 

INTER  interpolates  between  stored  A-weighted  overall  and  octave 
band  reference  levels  to  find  reference  intensities  in  fre- 
quency band  numbers  IQ  as  a function  of  the  average  speed  of 
vehicles  in  a group. 

FUNCTION  ACOS  (X) 

ACOS  is  a function  subroutine  which  calculates  the  arc  cos  of 
X between  0 and  tt. 

FUNCTION  AMAG  (XIV,  X2V) 

AMAG  is  a function  subroutine  which  calculates  the  distance 
between  the  points  XIV  and  X2V. 
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FUNCTION  ANGLE  (XIV,  X2V , X3V) 

ANGLE  is  a function  subroutine  which  calculates  the  angle  at 
point  X3V  that  is  subtended  by  the  line  between  XIV  and  X2V . 

The  angle  is  positive  with  values  between  0 and  tt. 

FUNCTION  BARFAC  (KF,  DELP) 

BARFAC  is  a function  subroutine  which  calculates  the  attenuation 
factor  due  to  barrier  diffraction  as  a function  of  the  frequency 
band  number  KF  and  of  the  path  length  difference  DELP  of  the 
diffracted  ray  path  relative  to  the  direct  ray  path. 

SUBROUTINE  BLOCN  (XIV,  X2V , X3V,  X4V,  X5V , X6V,  LOC) 

BLOCN  calculates  the  location  LOC  of  a barrier  or  a ground 
strip,  defined  by  the  end  points  X4V  and  X5V,  relative  to  the 
triangle  formed  by  two  points  XIV  and  X2V  on  a road  segment 
and  the  receiver  point  X3V  (see  Fig.  12).  The  point  X6V  is 
identical  with  X2V  if  both  points  X4V  and  X5V  are  outside 
the  triangle.  If  one  of  the  points  X4V  or  X5V  is  inside  the 
triangle,  X6V  is  the  intercept  of  a line  through  X3V  and  X4V 
with  the  road  segment  XIV,  X2V.  If  both  points  X4V  and  X5V 
are  inside  the  triangle,  then  the  point  X6V  is  the  intercept 
of  the  road  segment  XIV,  X2V  with  the  line  through  X3V  and 
either  X4V  or  X5V  depending  on  which  intercept  is  nearer  to 
the  point  XIV. 

LOC  = 0 denotes  that  the  line  X4V,  X5V  is  outside  the  tri- 
angle XIV,  X2V , X3V. 

LOC  = 1 denotes  that  the  line  X4V,  X5V  intersects  the  line 
XIV,  X3V. 

LOC  = 2 denotes  that  the  line  X4V,  X5V  intersects  the  line 
X2V,  X3V. 


8 7 


X5V 


XIV  X6V  X2V 


FIG. 12  NOTATION  FOR  SUBROUTINE  BLOCN 
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LOC  = 3 denotes  that  the  line  X4V,  X5V  intersects  both  lines 
XIV,  X3V  and  X2V , X3V. 

LOC  = 4 denotes  that  the  line  X4V,  X5V  is  completely  inside 
the  triangle  XIV,  X2V , X3V. 

FUNCTION  DEL  (XIV,  X2V , X3V,  X4V , HDIFF,  DN1) 

DEL  is  a function  subroutine  which  calculates  the  path  length 
difference  between  the  ray  emanating  from  the  source  point 
XIV  and  diffracted  over  the  barrier  X3 V,  X4V  towards  the 
receiver  X2V  and  the  direct  ray  from  XIV  to  X2V.  HDIFF  is 
the  height  of  the  direct  ray  above  the  barrier,  and  DN1  is 
the  path  length  traveled  by  the  direct  ray. 

FUNCTION  DSQR  (XIV,  X2V) 

DSQR  is  a function  subroutine  which  calculates  the  squared  length 
between  two  points  XIV  and  X2V. 

SUBROUTINE  ENDPT  (XIV,  X2V,  X3V,  X4V , X5V,  X6V,  KTRIG , IERR) 

ENDPT  calculates  the  new  end  point  of  a road  segment  from 
which  sound  is  uniformly  attenuated  during  propagation  towards 
the  receiver  X3V.  The  original  road  segment  has  the  initial 
point  XIV  and  the  end  point  X2V.  A barrier,  reflector,  or 
ground  strip  has  the  end  points  X4V,  X5V.  X6V  is  the  break 
point  on  the  road  between  segments  that  are  affected  and 
unaffected  by  the  barrier,  reflector,  or  ground  strip.  If 
the  entire  road  segment  XIV,  X2V  is  unaffected,  KTRIG  = 0 
and  X6V  is  identical  with  X2V . If  the  entire  road  segment 
XIV,  X2V  is  affected,  KTRIG  = 1 and  X6V  is  identical  with 
X2V.  If  a break  point  X6V  is  at  a distance  of  less  than 
0.51  ft  from  the  initial  point  XIV,  this  break  point  is  dis- 
regarded. If  the  initial  road  segment  XIV,  X6V  is  affected 
by  a barrier,  reflector,  or  ground  strip,  KTRIG  = 1,  the 
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break  point  X6V  is  moved  towards  XIV  by  0.5  ft,  and  the  new  end 
point  is  denoted  as  X6V.  If  the  initial  road  segment  is  un- 
affected, KTRIG  = 0,  the  break  point  X6V  is  moved  towards  XIV  by 
0.5  ft,  and  the  new  end  point  is  denoted  as  X2V. 

The  error  indicator  IERR  is  set  in  the  SUBROUTINE  MOVE. 

FUNCTION  GAMMA  (X) 

GAMMA  is  a function  subroutine  which  calculates  the  Gamma 
function  of  a real  constant  X. 

FUNCTION  GAUSS  (X) 

GAUSS  is  a function  subroutine  which  calculates  the  Gaussian 
probability  for  real  argument  X. 

FUNCTION  HEIGHT  (XIV,  X2V,  X3V,  X4V ) 

HEIGHT  is  a function  subroutine  which  calculates  the  height 
of  a ray  from  point  XIV  to  point  X2V  above  the  line  with  the 
end  points  X3V,  X4V. 

FUNCTION  I E PS  (XIV,  X2V,  X3V , X4V , DELI) 

IEPS  is  a function  subroutine  which  calculates  whether  the 
path  length  difference  DELI  is  sufficiently  similar  to  the 
path  length  difference  between  the  ray  emanating  from  XIV 
and  diffracted  over  the  barrier  X3V,  X4V  towards  the  receiver 
X2V  and  the  direct  ray  from  XIV  to  X2V.  If  the  two  path 
length  differences  are  sufficiently  similar,  IEPS  = 1;  other- 
wise, IEPS  = 0. 

SUBROUTINE  IMAGE  (XIV,  X2V,  X3V , X4V ) 

IMAGE  calculates  the  location  of  an  image  receiver  X4V  relative 
to  a receiver  X2V  and  a reflector  with  the  end  points  XIV, 

X2V . 
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SUBROUTINE  INTCPT  (XIV,  X2V,  X3V,  X4V , X5V ) 

INTCPT  calculates  in  the  x-y  plane  the  intercept  X5V  of  two  lines 
specified  by  the  two  pairs  of  points  XIV,  X2V  and  X3V,  X4V. 

FUNCTION  KCUT  (XIV,  X2V,  X3V,  X4V ) 

KCUT  is  a function  subroutine  which  checks  whether  two  line  seg- 
ments with  the  end  points  XIV,  X2V  and  X3V,  X4V,  respectively, 
cross  in  the  x-y  plane.  If  the  two  line  segments  cross, 

KCUT  = 1;  otherwise,  KCUT  = 0. 

FUNCTION  KPOS  (XIV,  X2V , X3V ) 

KPOS  is  a function  subroutine  which  checks  whether  the  point  X3V 
lies  between  two  points  XIV  and  X2V  on  a line.  If  it  does, 

KPOS  = 1;  otherwise,  KPOS  = 0. 

SUBROUTINE  MIDP  (XIV,  X2V,  X3V ) 

MIDP  calculates  the  center  point  X3V  between  two  points  XIV 
and  X2V. 

SUBROUTINE  MOVE  (XIV,  X2V,  X3V , DELTA,  IERR) 

MOVE  calculates  the  point  X2V,  on  the  line  specified  by  the  points 
XIV  and  X2V,  that  is  at  a distance  of  DELTA  feet  from  the  point 
XIV  in  the  direction  away  from  X3V.  MOVE  gives  the  error  message 
IERR  = 4 if  the  two  points  XIV  and  X3V  are  identical. 

SUBROUTINE  NRPT  (XIV,  X2V,  X3V  , X4V , DIST) 

NRPT  calculates  the  point  X4V,  on  a line  specified  by  the  points 
XIV,  X2V,  that  is  nearest  to  the  point  X3V.  It  also  calculates 
the  distance  DIST  between  the  points  X3V  and  X4V.  If  the  dis- 
tance is  zero,  DIST  = 1. 


91 


SUBROUTINE  NR1  (XIV,  X2V,  X 3V , X4V , DIST,  X5V  , DN1) 

NR1  calculates  the  point  X5V,  on  a line  segment  between  the  end 
points  XIV,  X2V,  that  is  nearest  to  the  point  X3V.  It  also  cal- 
culates the  distance  DN1  between  the  points  X3V  and  X5V.  The 
subroutine  employs  as  input  parameters  the  distance  DIST  between 
the  receiver  and  the  nearest  point  X4V  on  the  entire  line. 

SUBROUTINE  REPLCE  (XIV,  X2V) 

REPLCE  sets  the  vector  X2V  equal  to  the  vector  XIV. 

SUBROUTINE  SECTN  (XIV,  X2V,  X3V , X4V , X5V , X6V,  X 7 V ) 

SECTN  calculates  the  segment  X6V,  X7V  of  the  line  between  X4V  and 
X5V  that  is  bounded  by  the  two  lines  from  X3V  to  XIV  and  from 
X3V  to  X2V. 

FUNCTION  TAN  (X) 

TAN  is  a function  subroutine  which  calculates  the  tangent  of  X. 
SUBROUTINE  TRI  (XIV,  X2V,  X 3 V , X4V , X5V , KTRI ) 

TRI  determines  whether  or  not  the  point  X4V  lies  in  the  triangle 
formed  by  the  points  XIV,  X2V,  X3V.  If  X4v  is  in  the  triangle, 
KTRI  = 1;  otherwise,  KTRI  = 0.  TRI  also  calculates  the  inter- 
cept X5V  of  two  lines  specified  by  the  points  XIV,  X2V  and  X3V, 
X4v,  respectively. 

FUNCTION  XINGAM  (X,  TM ) 

XINGAM  is  a function  subroutine  which  calculates  the  incomplete 
Gamma  function  1-P  (TM,  X/2),  as  defined  in  Eq.  (39)  of  Sec.  2.4, 
by  making  use  of  approximations  with  the  Gaussian  probability 
function  as  given  in  Ref.  4.  If  the  argument  TM  is  too  small  for 
the  direct  application  of  an  approximation,  a recurrence  relation 
is  employed  that  increases  the  magnitude  of  TM. 
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FUNCTION  X P ROB  (A2,  AA3 , AA4 , A5 , X) 


XPROB  is  a function  subroutine  which  calculates  the  probability 
defined  in  Eq.  (39).  Three  parameters  correspond  to  the  notation 
in  Eq.  (39)  as  follows:  A2  = 3'/3",  AA3  = q'+q",  AA4  = q"  . The 

fourth  parameter  3’  and  the  normalized  intensity  I enter  into  A5 
and  X as  specified  in  the  main  program. 

FUNCTION  ZCOR  (XIV,  X2V,  X3V ) 

ZCOR  is  a function  subroutine  which  calculates  the  z-coordinate 
of  the  point  X3V  so  that  the  point  X3V  lies  on  the  line  speci- 
fied by  the  points  XIV  and  X2V. 
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3.6  Computer  Program  Listing 


C MAIN  PROGRAM  12/14/71 

DIMENSION  SIG(2) ,XKAP(4) 

DIMENSION  XR 1 ( 3 ) ,XR2 ( 3) 

COMMON  /INOU/INPT, TOUT 
COMMON  /BLK2/NQ 
COMMON  /INPT 1 /RDIN ( 15) 

COMMON  /INPT2/NR,NB,NG 

COMMON  /INPT3/XRC(15),YRC(15),ZRC(15),NRC 
COMMON  /DRIV2/NQS( 10,2) ,NF 
COMMON  /DRIV3/CQ(2) ,XLE(9) 

COMMON  /STORE4/XMPH(  10,5,2),  VEXPH ( 1 0 , 5 , 2 ) 

COMMON  /STORE7/  XLB ( 5 ) , STB ( 1 5 ) 

COMMON  /DRIV4/XKA(5) 

COMMON  /ST0RE3/RX( 10, 1 1 ) ,RY ( 10, 1 1 ) ,RZ ( 10, 1 1 ) ,NRSM1  ( 10) 
COMMON  /GE01/IBAR,XSEG,IGRA 
DATA  LZ/2HZZ/ 

INTEGER  TITLE( 15) 

EQUIVALENCE  ( RDIN ( 3 ) , SIG ( 1 ) ) 

1 N Q = 2 
NG  = 0 
N B = 0 

CALL  FILES 

5 READ(INPT,  1005)  ( TITLE ( I ) ,1=1,  15) 

IF  ( TITLE ( 1 ) ,EQ.LZ)GO  TO  105 
WRITE(IOUT,2009) 

WRITE(IOUT, 1005)  ( TITLE ( I ) ,1=1, 15) 

CALL  INPUT 
NF=RDIN ( 2 ) 

WRITE(IOUT, 1002) (TITLE(I) ,1=1, 15) 

WRITE ( IOUT  f 1 003  ) 

DO  11  1=1, NQ 

CQ ( I ) =EXP ( .5*(SIG(I)*.23026)**2) 

11  CONTINUE 

DO  60  1=1, NBC 
DO  15  J=  1 , NF  ■ 

X L E ( J ) =0  * 

15  CONTINUE 

DO  16  J = 1 , 4 
XKA ( J )=0. 

16  CONTINUE 

DO  50  M = 1 , N R 
DO  20  I Q=  1 , NQ 

IF (NQS(M,IQ) , EQ , 0 ) G Q TO  20 
CALL  INTER ( M, IQ ) 

20  CONTINUE 

XR  1 ( 1 )=RX(M, 1 ) 

XR  1 (2)=RY(M,  1 ) 

X R 1 ( 3 ) =R Z ( M , 1 ) 

NLIM=NRSM 1 (M) 
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DO  40  N=2,NLIM 
X R 2 ( 1 ) =RX ( M, N ) 

X R 2 ( 2 ) =RY ( M,  N ) 

XR2(3)=RZ(M,N) 

CALL  GE0MRY(XRC(I),YRC(I),ZRC(I),XR1,XR2,IERR,M) 
IF(IERR,EQ, 1 ) G 0 TO  65 

IF ( IERR ,EQ. 3 ) GO  TO  66 
IF ( IERR  ,EQ,  4 ) GO  TO  60 
CALL  REPLCE(XR2»XR1 ) 

40  CONTINUE 
50  CONTINUE 

DO  55  J = 1 , N F 

XLE(J)=100,+10.*ALOG  |0(XLE(  J)  ) 

55  CONTINUE 

WRITE (I OUT, 1004)  I , XR  C ( I ) ,YRC(I) , Z R C ( I ) 

IF(NF.EQ.I)  GO  TO  56 
WRITE ( IOUT,  200  1 ) 

WRITE (IOUT, 2002)  (XLE(II) ,II=2,NF) 

56  DO  57  J = 1 , 4 
A J = J 

XKA( J)=XKA( J)+10.**(AJ/10.*(XLB(I)-10I0,+AJ*STB(I) **2/8,7) ) 

57  CONTINUE 

XLA=100.+10,0*ALOG10(XKA(1)) 

DO  59  J = 2 , 4 

XKAP( J)=XKA( J)/XKA(  1 )**J 
5D  CONTINUE 

SIGL  = 4.35*SQRT(AL0G(  1 .+XKAP(2)  ) ) 

XLNP=XLA+2,56*SIGL 
XL50=XLA-SIGL**2/8,7 
XL10=XL50+1 ,28*SIGL 
XL90=XL50-1 ,28*SIGL 

TEMP 1=XKAP( 4 )+3.*XKAP( 2) *XKAP( 3 ) 
TEMP2=3.*(XKAP(3)-2.*XKAP(2)**2) 

R=2.*TEMP2/TEMP1 

S=2.*XKAP(2)-XKAP(3) * TEMP 2/TEMPI 
IF ( 1 ,/R*S  ,GE. 1 . ) GO  TO  200 
TEMP  1 =SQRT ( 1 ,-R*S) 

A 1 = ( 1 .-TEMPI  )/R 
A2=2./S*A1-1 . 

AA3=R/S*(2,/R-XKAP(2) ) 

AA4=(AA3+(R-AA3) /TEMPI )/2. 

TEMP 1=GAMMA ( AA  3 + 1 . ) 

TEMP  1=A2**AA4/TEMP1 
XI=1  . 

FLAG=1  . 

ZI=XLA+1 . 

170  Z 1 =XI/ A 1 

TEMP2  = TEMP1*Z1**AA3/EXP(Z1  ) 

PRO=XPROB(A2,AA3,AA4,TEMP2,Z1) 

ZIsXLA-FLAG 
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PR=GAUSS ( (ZI-XLB(I) )/STB(I) ) 
IF(PRO.GT,PR)PRO*PR 
IF(FLAG.GT.0.0)GO  TO  180 
IF(PRO.GE,0#9)GO  TO  190 
XI*XI*1 .25867 
GO  TO  170 

180  IF(PRO.GT#0,9)XI  = XI/  1 .25867 
IF(PRO.LE.0.9)FLAG=- 1 . 

IF ( PRO. LE,0. 9 )XIsXI* 1.25867 
GO  TO  170 

190  IF (ZI.GT.XL50)XL10=ZI 


200  WRITE ( IOUT, 2003 ) 

WRITE (IOUT. 2004 )XLA#XLNP,XL90,XL50, XI 10 

IF ( XL B ( I ) » EQ  « 0 a 0 ) GO  TO  60 
DXL A=XL A-XLB ( I ) -STB(I)**2/8,7 
DXLNP  = XLNP  + DXLA--2.56*STB(I)f.XLA 
DXL90=XL90-XLB(I)+1 ,28*STB(I) 

DXL50=XL50-XLB(I) 

DXL10=XL10-XLB(I)-1 .28*STB(I) 

WRITE ( IOU T, 2005  )DXLA,DXLNP.DXL90,DXL50, DXL  10 

60  CONTINUE 
GO  TO  5 

65  WRlTEdOUT, 1006)1, M,N,XBAR, IS  EG 
GO  TO  5 

66  WRITE(I0UT,  1009  )1, M,N 
GO  TO  5 

70  WRITE(I0UT, 1008)1, M,N,IGRA 
GO  TO  5 

105  CALL  RELEAS(INPT) 

CALL  RELEAS ( TOUT ) 

GO  TO  1 

1002  FORMAT(//15A2) 

1003  FORMAT ( 1 1H0PECEIVER  # 1 0 X , 3 HXR C 9 X , 3H X RC 9X , 3H Z RC ) 

1004  F0RMAT(/4X,I3,5X,3F12. 1 ) 

1005  F0RMAT(15A2) 

1006  FORMAT ( 44H0ILLEGAL  BARRIER  INTERSECTS  RQADWAY  FOR  REC#I2 , 2X , 2HR#I2 
1,3HRS#I2,2X,2HB#I2,2X,3HBS#I2) 

1008  FORMAT { 49H0ILLEGAL  GROUND  STRIP  INTERSECTS  ROADWAY  FOR  REC#I2»2X,2 
1HR#I2,3HRS#I2,2X,4HAGS#I2) 

1009  FORMAT ( 26H0TOO  MANY  REFLECTIONS , RCV# , 12 , 4H  R#,I2,4H  S#,I2) 

2001  FORMAT(/14X,2H635X,3H1254X,3H2504X,3H5003X,4H10003X,4H20003X, 

X 4H40003X, 4H8000  ) 

2002  FORMAT( 10X,8F7.  1 ) 

2003  FORMAT(/10X,5HLE(A),5X,3HLNP,5X,3HL90,5X,3HL50,5X,3HL10) 

2004  FORMAT ( 8X , 5F8  , 1 ) 

2005  F0RMAT(3X,5HINCR.  ,5F8.  1 ) 
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o n 


2009  FORMAT ( 1H 1 , 15X, 24HTRAFFIC  NOISE  PREDICTION//) 

END 

$1 

SUBROUTINE  FILES 
COMMON/INOU/INPT, I OUT 
IN  PT= 1 
IOUT=23 
RETURN 
END 
SI 

SUBROUTINE  INPUT 

COMMON  /DRIV2/NQS ( 10,2  ) ,NF 

INTEGER  ALPH A ( 15) 

DIMENSION  VTEMP ( 5 ) 

COMMON  /STORE1/BX(10,5),BY(10,5),BZ(10,5) , I BLAST (10),NBSM1(10) 
COMMON  /STORE2/XXG1  (10,2),YYG1(10,2),ZZG1(10,2),BGS(10),IDUM(10) 
COMMON  /STORE3/RX(10,11),RY(10,11),RZ(10,11),NRSM1(10) 

COMMON  /STORE4/XMPH ( 10,5,2)  , VEXPH ( 1 0 , 5 , 2 ) 

COMMON  /STORE7/  XLB ( 1 5 ) , STB ( 1 5 ) 

COMMON  /IN0U/INPT,I0UT 
COMMON  /BLK2/NQ 
COMMON  /INPT  1 /RDIN (15) 

COMMON  /INPT2/NR,NB,NG 

COMMON  /INPT3/XRC(  15)  ,YRC( 15) ,ZRC(  15) ,NRC 
EQUIVALENCE  ( B DIN ( 1 ) , XNIGHT ) 

DATA  IA/2HA  / 

DATA  IG/2HG  / 

DATA  IR/2HR  / 

DATA  IT/2HT  / 

DATA  LAST/2HL  / 

5 READ(INPT,  1000) IGO,1 1 ,12 

GO  TO(  10,20,30,36,40,50  ) ,IGO 
C GARBAGE  DATA. .. .PROGRAM  INITIALIZATION  PARAMETERS 

10  WRITE(IOUT,2000) 

1 1 READ (INPT, 1001 ) VALUE, IDN,IL AST, (ALPHA (I) ,1*1, 15) 

RDIN ( IDN ) =VALUE 

WRITE(IOUT,2001)VALUE,IDN,(ALPHA(I),I=1,15) 

IF(lLAST.NE.LAST)GO  TO  11 
GO  TO  5 
C 


C 


VEHICLE  DATA 

20  NR=I1 

DO  28  J=  1 , NR 
NSEC= 1 

DO  21  K=1,NQ 

21  NQS ( J , K ) °0 

22  READ(INPT,  1002)  VEH , XMH , ITY , ILAST 
NQS( J,ITY)=NQS( J,ITY)+1 

NQC=NQS ( J , ITY ) 

VEXPH (J,NQC, ITY) =VEH/XMH/5280  . 
XMPH( J,NQC,ITY)*XMH 
IF (ILAST, NE.LAST)GO  TO  22 
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C ROADWAY  DATA  SECTIONS 

24  READ(INPT,1003)RX(J,NSEC),RY(J,NSEC),RZ(J,NSEC),ILAST 
IF (ILAST.EQ,LAST)GO  TO  25 

NSEC=NSEC+1 
GO  TO  24 

25  IF ( N S EC - 1 , N E . 0 ) G 0 TO  26 
WRlTE(IOUT,2010) 

CALL  EXIT 

26  NRSM  I ( J)=NSEC-1 
28  CONTINUE 

GO  TO  5 

C BARRIER  DATA  SECTIONS 

30  N B = I 1 

IF(NB.EQ.0)GO  TO  5 
DO  35  J=1,NB 
NSFC= 1 

31  READ(INPT,1003)8X(J, NSEC), BY(J, NSEC ),BZ(J, NSEC), IBLAST(J) 
IF(IBLAST( J)  .EQ.IA.OR,IBLAST( J)  ,EQ,IR)GO  TO  32 
NSEC=NSEC+ 1 

GO  TO  31 

32  IF(NSEC-1 ,NE.0)GO  TO  33 
WRITE ( IOUT, 20 1 1 ) 

CALL  EXIT 

33  NBSM1(J)=N SEC-1 

35  CONTINUE 
GO  TO  5 

C ABSORBING  GROUND  STRIPS 

36  NG  = I 1 

IF(NG.EQ.0)GO  TO  5 
DO  37  1=1, NG 

READ(INPT,1004)XXG1(I,1),YYG1(I,1),ZZG1(I,1),BGS(I) 

READ(INPT,1003)XXG1(I,2),YYG1(I,2),ZZG1(I,2),IDUM(I) 

IF ( IDUM ( I ) ,EQ.IG)IDUM(I)=1 
IF ( IDUM ( I ) ,EQ,IT)IDUM(I)=2 

37  CONTINUE 
GO  TO  5 

C RECEIVER  DATA 

40  NRC=I1 

DO  41  1=1, NRC 

READ(INPT,1005)XRC(I),YRC(I),ZRC(I) ,XLB(I) ,STB(I) 

ZRC(I)=ZRC(I)+XNIGHT 

41  CONTINUE 
GO  TO  5 

C 

50  DO  65  J = 1 , N R 

WRITE(IOUT,2002) J 
DO  55  K = 1 , N Q 

IF(NQS( J,K) ,EQ.0)GO  TO  55 
NLIM=NQS( J,K) 

DO  54  1=  1 , NLIM 
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5U  VTEMP(I)=VEXPH( J,I,K) *XMPH( J,I,K)*5280, 

WHITE  (IOUT,  2004  )K, (I,VTEMP(X)  , XMPH ( J , I , K ) , I* 1 , NQS ( J , K ) ) 

55  CONTINUE 

WRITE(IOUT,2013) 

WRITE(IOUT,2005)RX(J,1),RY(J,1),RZ(J,1) 

NSEC=NRSM1 ( J)+1 
DO  60  1 = 2 , NSEC 

WRITE (IOUT, 2006)1, RX(J, I), RY(J, I), RZ(J,I) 

60  CONTINUE 
65  CONTINUE 

IF ( NB , EQ , 0 ) GO  TO  80 
DO  75  J=1,NB 

WRITE ( IOUT ,2007) J,IBLAST( J) 

WRITE(IOUT,2005)BX(J,1),BY(J,1),BZ(J,1) 

NSEC=NBSM 1 ( J ) 4 1 
DO  70  1=2 , NSEC 

WRITE (IOUT, 2006 ) I,BX ( J,I) ,BY< J,I) ,BZ(J,I) 

70  CONTINUE 
75  CONTINUE 

80  IF(NG.EQ.0)  GO  TO  90 
DO  85  1=1, NG 
I F ( I D UM ( I ) .EQ. 1 )IDM  = IG 
IFflDUM(I) , EQ , 2 ) IDM=IT 

WRITE (IOUT, 2012)1, IDM, XXG 1(1, 1 ) ,YYG  1(1, 1),ZZG1(I,1), 

1 BGS(I),XXG1(I,2),YYG!(I,2),ZZG1(I,2) 

85  CONTINUE 
90  WRITE ( IOUT. 2008 ) 

DO  95  1=1, NRC 

WRITE (IOUT, 2014 )I,XRC (I) , YRC(I) ,ZRC(I) ,XLB(I) , STB ( I ) 

95  CONTINUE 
RETURN 

1000  F0RMAT(3I5) 

1001  FORMAT(E10.0,I5,4X,A1,10X,15A2) 

1002  FORMAT(2E10.0,I5,5X,A'i  ) 

1003  FORMAT(3E10.0, A1 ) 

1004  FORMAT(4E10,0) 

1005  FORMAT(3E10.0,2F10. 1 ) 

2000  FORMAT ( 34H0PROGRAM  INITIALIZATION  PARAMETERS) 

2001  FORMAT(1X,E12.5,I10,5X,15A2) 

2002  FORMAT ( 10H0ROADWAY  #,I3) 

2004  FORMAT ( 10H0NUMBER  OF  1 3X, 5HVEH/H8X, 3HMPH/5H  TYPE, 12, 4H  VEH/(3x,I2 
1,  1 5X , 2E  1 3 , 4 ) ) 

2005  FORMAT(7H0NUMBER,5X,  1 H X 1 2 X , 1HY12X, 1HZ/4X, 1H1,2X,3E13.4) 

2006  F0RMAT(3X,I2,2X,3E13.4) 

2007  FORMAT ( 10H0BARRIER  #13 , 2X , 1 H ( A 1 , 1 H ) 1 0X , 1 9HBARRIER  COORD  IN  FT) 

2008  FORMAT(9H0RECEIVER14X,20HRECEIVER  COORD  IN  TT/7H  NUMBER5X , 1 HX  1 2X , 1 
1HY,  12X,  1HZ, 12X,5HLBACK,7X,7HSTDEVLB) 

2010  FORMAT (27H0INSUFFICIENT  ROAD  SECTIONS) 

2011  FORMAT (30H0INSUFFICIENT  BARRIER  SECTIONS) 

2012  FORMAT ( 18H0ABSORBING  STRIP  #13 , 2X , 1 H ( A 1 , 1 H ) //5H  PT  #7X,  1 HX  1 2X , 1HYi 
A2X, 1HZ12X,5HWIDTH/4X, 1H12X,4E13.4/4X,1H22X,3E13.4) 

2013  FORMAT(22X, 18HS0URCE  COORD  IN  FT) 

2014  FORMAT(3X,I2,2X,3E13.4,2F10, 1 ) 

END 

$0 
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C GEOMETRY  ROUTINE  12/13/71 


SUBROUTINE  GEOMR Y ( XR , YR , ZR , XR  1 0 , XR 2 0 , IERR , MR ) 

DIMENSION  XG1(3),XG2(3),XG3(3),XG4(3) 

DIMENSION  B1(3,2,10),R1(3,2,10),RB1(3,2,10),TA1(3,2,10) 
DIMENSION  KBCODE ( 10) ,KNUMB(  10) ,KRNUMB(  10) , KRDNUM ( 10) 

DIMENSION  KGCODE(10),BGT(10),IKIN(10) 

DIMENSION 

D DELP0 ( 2 ) , DELP 1 ( 2 ) , DELP2 ( 2 ) , FB ( 9 , 2 ) , PELR ( 2 , 1 0 ) , FG ( 9 ) , HGA ( 2 ) , 

X XB1 (3) ,XB2(3) ,XDB1 (3) , XDB2 ( 3 ) , XDB3 ( 3 ) ,XDB4(3), 

X XRB 1 (3) , XRB2 ( 3) , XRB3 ( 3 ) , XRB4 ( 3 ) , XRC ( 3 ) , XRCI ( 3 ) , 

X XR 1 (3) ,XR2(3) ,XR1I(3) ,XR2I(3) . XR  1 0 ( 3 ) , XR2  0 ( 3 ) , 

X X K ( 3 ) ,XKI(3) ,XJ(3) ,XNPT(3) *XNPTI(3) ,XNPTJ(3)  , 

X XN1(3),XN2{3) ,XN1I(3) ,XIMG(3, 1 0 ) , XR2D ( 3 ) , XR2G ( 3 ) , 

X XL A ( 5 , 2 ) , 

Z Z S ( 2 ) 

COMMON /INOU/INPT, TOUT 
COMMON/BLK2/NQ 
COMMON/INPT1/RDIN ( 15 ) 

C0MM0N/INPT2/NR,NB#NG 
COMMON/DRI V2/NQS ( 10,2) , NF 
COMMON/DRI V 3 /CQ(2) ,XLE(9) 

COMMON  /STORE4/XMPH ( 1 0 , 5 , 2 ) , V EX  PH ( 1 0 , 5 , 2 ) 

COMMON/STORE1/BX(10,5),BY(10,5$,BZ(10.5),IBLAST(10),NBSM1(10) 
COMMON/STORE2/XXG1 (10,2),YYG1(10,2),ZZG1(10,2),BGS(10),IDUM(10) 
C0MM0N/DRIV4/XKA ( 5 ) 

COMMON/GEO  1 /IB AR,ISEG,IGR A 
COMMON/INTER1/XLREF( 10,5,9,2) 
equivalence  ( RDIN ( 5 ) , ZS ( 1)  ) 

DATA  IA/2HA  / 

DATA  IR/2HR  / 

DATA  HGA(  1) , HG A ( 2 ) / 1 0 . , 30  . / 

XRC ( 1 ) =XR 
X R C ( 2 ) = Y R 
XRC ( 3 ) =ZR 
IERR=0 

CALL  NRPT(XR10,XR20,XRC,XNPT»DIST) 

ANG1=ANGLE(XR10,XNPT,XRC) 

ICODE= 1 

WRITE ( IOU T,  1000  )ICODE,XR  10, XR20,XRC,XNPT 
PRELIMINARY  SELECTION  OF  BARRIERS 
NDIFF=0 

IF ( N B # EQ  „ 0 ) GO  TO  16 
KBAR=0 

DO  15  IBAR=1,NB 
KAR  = IBLAST(IBAR  ) 

XB1 ( 1 )=BX(IBAF,  1 ) 

XB1 ( 2 ) =B Y ( I8AR , 1) 

X B 1 ( 3 ) =B Z ( IB A R , 1 ) 

NLIM=NBSM 1 (IBAR ) +1 
DO  15  ISEG=2 , NLIM 
X B 2 ( 1 )=BX(IBAR,ISEG) 

XB2(2)=BY(IBAR,ISEG) 

XB2(3)=BZ(IBAR,ISEG) 

KBAR=KBAR+1 
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C IC0DE=2 

C WRITE(IOUT,1001)ICODE,KBAR,XB1,XB2 

IF(KCUT(XR10#XR20»XB1#XB2) .NE.IJGO  TO  10 

IERR= 1 

RETURN 

10  CALL  BLOCN(XR10,XR20,XRC,XB1,XB2,XK,LOC) 
IF(LOC.EQ.0)GO  TO  11 

ndiff=ndiff+i 

KNUMB(NDIFF)=KBAR 

KBCODE ( NDIFF ) =LOC 

CALL  REPLCE(XB1,B1 ( 1, 1, NDIFF)  ) 

CALL  REPLCE(XB2,B1 ( 1,2, NDIFF)  ) 

11  CALL  REPLCE ( XB2  # XB  1 ) 

15  CONTINUE 

C PRELIMINARY  SELECTION  OF  STRIPS 

16  KG A=0 

IF(NG.EQ.0)GO  TO  20 
DO  18  IGRA*1,NG 
XG1 (1 ) =XXG 1 (IGRA,  1 ) 

XG1 ( 2 ) =Y YG 1 (IGRA,  1 ) 

XG1 ( 3 ) = Z Z G 1 (IGRA,  1 ) 

X G 2 ( 1 )=XXG1 ( IGR A , 2 ) 

XG2 ( 2 ) =YYG 1 ( IGR A , 2 ) 

XG2(3)=ZZG1 ( IGR A , 2 ) 

ICODE=3 

WRITE(IOUT,1001)ICODE,IGRA,XG1,XG2 
IF(KCUT(XR10,XR20#XG  i,XG2) ,NE. 1 )GO  TO  17 
IERR=2 
RETURN 

17  CALL  BLOCN(XR10,XR20,XRC,XG1,XG2,XK,LOC) 
IF (LOC.EQ.0)GO  TO  18 

kga=kga+i 

KGCODE(KGA)=LOC 
CALL  REPLCE(XG1,TA1  ( 1#  1 # KG A ) ) 

CALL  REPLCE(XG2,TA1  ( 1 , 2 * KG A ) ) 

BGT(KGA )=BGS(IGRA) 

IKIN (KGA)=IDUM(IGRA) 

18  CONTINUE 

C DIFFRACTION  OF  DIRECT  RAY 
20  CALL  REPLCE(XR10,XR1 ) 

25  CALL  REPLCE(XR20,XR2) 

DO  30  IQ=1,NQ 
DELP0 ( IQ ) 3- , 2 
30  CONTINUE 
C ICODE=U 

C WRlTE(IOUT, 1000)ICODE,XR1,XR2 

IF(NDIFF.EQ.0)GO  TO  71 
ITRIG=0 

DO  70  KDIFF=1 , NDIFF 
KBAP.  = KNUMB(KDIFF) 

KCD=KBC0DE(KDIFF) 

CALL  REPLCE(B1 ( 1, 1,KDIFF) ,XB1 ) 

CALL  REPLCE(B1 ( 1,2,KDIFF),XB2) 
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C IC0DE=5 

C WRITE(IOUT,1001)ICODE,KBAR.XR1,XR2,XB1#XB2 

IF ( KCD .EQ. 3 ) GO  TO  40 

CALL  ENDPT(XR1,XR2,XRC,XB1,XB2,XK,KTRIG,IERR) 

IF(IERR.EQ.U)  RETURN 
IF(KTRIG.EQ.0)GO  TO  70 
GO  TO  45 

40  CALL  REPLCE(XR2,XK) 

45  MODD=0 
C ICODE=6 

C WRITE(IOUT,1001)ICODE,KCD,XR2,XK 

DO  60  IQ=1,NQ 

CALL  NR1(XR1»XK,XRC,XNPT,DIST,XN1,DN1) 

ZN  10  = XN  1(3) 

X N 1 ( 3 ) =ZN  1 0 + ZS ( IQ ) 

HDIFF=H EIGHT ( XN1 , XRC, XB 1 , XB2 ) 

IF ( IQ. ME. 1 ) GO  TO  50 
IF (HDIFF  . G T , 2 0 . )GO  TO  70 
50  DELP=DEL(XN1,XRC,XB1,XB2,HDIFF,DN1) 

IF(DELP.GT,DELP0(IQ) )GO  TO  53 

52  IF(MODD.EQ, 1 )GO  TO  65 
GO  TO  70 

53  IF  (IQ.NE.1)  GO  TO  59 
DR1  = AMAG ( XRC , XR1 ) 

IF  (ABS(DRI-DNI) . LT  • 1. ) GO  TO  54 
HDIFA  = HEIGHT  ( XR  1 , XRC , XB  1 , XB2 ) 

DELPA  = DEL  ( XR  1 , XRC , XB  1 , XB2 , HDIFA , DR  1 ) 

D E L M = (DELPA  + DELP)/2. 

C ICODE=107 

C WRITE ( TOUT, 1000) ICODE,DELP, DELPA 

IF ( ( ABS (DELPA-DELP ) -0. 1-DELM/50. * ( 1 ,+DELM) ) ,LE.0. ) GO  TO  55 
CALL  MIDP  (XR1,XN1,XK) 

IC0DE=7 

WRITE(I0UT,  1000) IC ODE, XRl,XN1,XK 
DEEP  = DELPA 

54  IF  (IEPS(XK,XRC,XB1,XB2,DELP) ,EQ.0)  GO  TO  58 
CALL  MIDP  (XR1,XK,XK) 

C ICODE=8 

C WRITE (IOUT,  1000)ICODE,XR1 ,XK 

GO  TO  54 

55  DRK=  AMAG  ( X R C , X K ) 

IF  (ABS(DRK-DN1  ) .LT.  .)  GO  TO  58 

56  IF ( IEPS  (XK,XRC,XB1,XB2,DELP)  .EQ.0)  GO  TO  58 
CALL  MIDP(XN1,XK,XK) 

GO  TO  56 

58  IF  (DELP,LE,DELP0(  1 ) ) GO  TO«52 

59  DELP0  (IQ)  = DELP 

C ICODF=  1 0 
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C WRITE(IOUT,1000)ICODE,DELP0(1),DELP0(2) 

MODD= 1 
60  CONTINUE 

IF(DELP0(2) ,LT. 1 2 • 5 ) GO  TO  65 
I TR IG= 1 

IF(KCD.NE.3,OR.KCD.ME.2)GO  TO  65 
CALL  REPLCE(XR20,XR2) 

GO  TO  66 

65  CALL  REPLCE(XK,XR2) 

66  IF (ITRIG.EQ, 1 ) GO  TO  71 

70  CONTINUE 

C ICODE=1  1 

C WRlTE(IOUT, 1000)ICODE,XR2 

71  ANG=ANGLE ( XR 1 , XR2,XRC) 

DELP  = DELP0 ( 1 ) 

FB( 1, 1 )=BARFAC( 1,DELP) 

C ICODE= 1 2 

C WRiTE(IOUT,  1002)ICODE,FB(  1,  1 ) 

CALL  REPLCE ( XR2, XR2D ) 

C GROUND  ABSORPTION 

72  DO  73  KF=1,NF 
FG ( KF ) = 1 . 

73  CONTINUE 

C ICODE=  1 3 

C WRITE(IOUT#1000)ICODE,XR1,XR2,XR2D 

IF(DELP0(  1 ) ,GT.-0.2)GO  TO  78 
IF(KGA.EQ.0)GO  TO  78 
DO  77  IGRA=1,KGA 
LOC=KGCODE(IGPA) 

CALL  REPLCE (TA1 ( 1 , 1 ,IGRA  ) , XG1  ) 

CALL  REPLCE(TA1(1,2#IGRA),XG2) 

BG  = BGT ( IGR A ) 

C IC0DE=14 

C WRITE ( IOUT, 1001 )ICODE,LOC,XG1,XG2 

IF(LOC.EQ.3)GO  TO  74 

CALL  ENDPT(XR1#XR2,XRC#XG1#XG2,XK#KTRIG,IERR) 
IF ( IERR . EQ  » 4 ) RETURN 
IF (KTRIG.EQ,0)GO  TO  77 
GO  TO  75 

74  CALL  REPLCE(XR2,XK) 

75  CALL  NR1(XR1,XK,XRC*XNPT,DIST,XN1#DN1) 
HDIFF=HEIGHT(XN1,XRC,XG1#XG2) 
IF(HDIFF.GT.HGA(IKIND) )GO  TO  77 

CALL  SECTN(XR1,XK,XRC,XG1,XG2,XG3,XG4) 

DL  = 1 ,57/(  1 ,/BG+l ./AMAG(XG3,XG4)  ) 

DO  76  IK= 1 # NF 
P P = IK 

IF(IK.EQ. 1 )PP=5 

IF(IKIND.EQ,1)A=(.00i6*PP-0.0028)*DL 

IF( IKIND.EQ.2) A=2. ** (PP/3 . ) / 1 3 1 0 . *DL 
IF( A.GT.3, ) A=3. 

FG(IK)=FG(IK)/10.**A 
IF(FG(IK).LT.1.E-3)FG(IK)=1.E-3 
CALL  REPLCE ( XK  # XR2 ) 

76  CONTINUE 
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on  on  on 


C IC0DE=15 

C WRITE(IOUT,1000)ICODE,FG(1),FG(9) 

77  CONTINUE 

78  CALL  REPiCE(XB2,XR2G) 

C PRELIMINARY  SELECTION  OF  REFLECTORS 
NREF=0 

IF ( NB  , EQ , 0 ) GO  TO  91 

N R F D F = 0 

KBAR=0 

DO  85  IBAR= 1 , NB 
KAR=IBLAST(IBAR) 

XRB1  ( 1 )=BX(IBAR,  1 ) 

XRR1  (2)=BY(IBAR# 1) 

XRB1  ( 31=BZ(IBAR, 1 ) 

NLIM  = NBSM1  (IBAR)-M 
DO  85  ISEG=2,NLIM 
X R B 2 ( 1 )=BX(IBAR,ISEG) 

XRB2(2)=BY(IBAR,ISEG) 

XRB2(3)=BZ(IBAR,ISEG) 

kb ar=kbar+ i 

C IC0DE=  1 6 

C WRITE(IOUT,1001)ICODE,KBAR,XR1,XR2,XRB1,XRB2 

CALL  IMAGE(XRB1,XRB2,XRC,XRCI) 

CALL  NRPT(XR1,XR2,XRCI,XNPTI,DISTI) 

ICODE= 1 7 

WRITE ( IOUT,  1000 )ICODE,XRCI,XNPTI 
ANGIMG=ANGLE(XRB1,XRB2,XRCI) 
FCTR=(DIST*ANGIMG)/DISTI/ANG/FB( 1,1) 

IF ( FCTR .LT, 0. 1 ) G 0 TO  80 
NRFDF=NRFDF+1 
KRDNUM(NRFDF)=K8AR 
CALL  REPLCE(XRB1,RBl  ( 1,  1,NRFDF)  ) 

CALL  REPLCE ( XRB2, RP 1 ( 1 , 2, NRFDF)  ) 

IC0DE=  1 8 

WRITE(IOUT,1001)ICODE, NRFDF, XRB1,XRB2 
IF(KAR.NE.IR)GO  TO  80 

CALL  BL0CN(XR1,XR2,XRCI,XRB1,XRB2,XK,L0C) 

IF ( LOC . EQ , 0 ) G 0 TO  80 

nref=nref+i 
KRNUMB(NREF)=KBAR 
CALL  REPLCE(XRB1*R1  ( l#  1 * N R E F ) ) 

CALL  REPLCE ( XRB2,R 1 ( i,2,NREF)  ) 

ICODE=  1 9 

WRITE (I OUT, 1001 )IC0DE,NREF,XRB1,XRB2 
80  CALL  REPLCE(XRB2,XRB'i  ) 

85  CONTINUE 
C BEGIN  REFLECTOR  PROBLEM 
91  IDXR=0 

IF(NRSF.EQ.0)GO  TO  180 
DO  175  KREFs 1 , NREF 
KBAR  1=KRNUMB ( KREF ) 

CALL  RF.PLCE(R1  ( 1,  1 , K R E F ) ,XRB1  ) 

CALL  REPLCE(R1 ( 1,2, KREF)  ,XRB2) 
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C ICODE=20 

C WRiTE(IOUT, 1001 )IC0DE,KBAR1,XRB1,XRB2,XR1,XR2 

CALL  IMAGE(XRB1,XRB2,XRC,XRCI) 

CALL  ENDPT(XR1,XR2,XRCI,XRB1,XRB2,XK,KTRIG,IERR) 
IF ( IERR ,EQ, 4 ) RETURN 
C IC0DE=21 

C WRITE(IOUT,1001)ICODE,KTRIG,XRCI,XK 

IF (KTRIG,EQ,0)GO  TO  175 
CALL  NRPT(XR1,XR2,XRCI,XNPTI,DISTI) 

CALL  NR1(XR1,XK,XRCl,XNPTI,DISTI,XN1I,DN1I) 

C ICODE=22 

C WRlTE(IOUT,1000)ICODE,XNPTI,XN1I 

XN1I(3)=XN1I(3)+ZS(2) 

HDIFF =HEIGHT(XN1I,XRCI,XRB1,XRB2) 

IF ( HDIFF .GT,-2.0)GO  TO  175 
CALL  FEPLCE(XK,XR2) 

CALL  SECTN(XR1,XR2,XRCI,XRB1,XRB2,XRB3,XRB4) 

C ICODE=23 

C WRITE(IOUT,1000)ICODE,XR1,XR2,XRB3,XRB4 

MDIFF=0 
DO  95  IQ=1,NQ 
DELP1  ( IQ ) =-0 . 2 
DELP2 (IQ) =-0,2 
95  CONTINUE 

C DIFFRACTION  BEFORE  REFLECTION 
IF ( NRFDF , EQ , 0 ) GO  TO  1 15 
DO  110  KRFDF=1 , NRFDF 
KBAP-2  = KRDNUM  ( KRFDF  ) 

CALL  REPLCE(RB1(1,1, KRFDF), XDB1) 

CALL  REPLCE(RB1 ( 1,2, KRFDF) ,XDB2) 

C ICODE=24 

C WRITE(IOUT,1001)ICODE,KBAR2,XDB1,XDB2 

IF(KBAR2.E0,KBAR1 )GO  TO  110 

CALL  BLOCN(XRB3,XRB4,XRCI,XDB1,XDB2,XK,LOC) 

IF (LOC.NE,0)GO  TO  1 10 

CALL  ENDPT(XR1,XR2,XRCI,XDB1,XDB2,XK,KTRIG,IERR) 
IF ( IERR . EQ , 4 ) RETURN 
C ICODE=25 

C WRiTE(IOUT, 1001 )ICODE,KTRIG,XK 

IF (KTRIG,EQ,0)GO  TO  1 10 

CALL  NR  1 (XR1,XK.XRCI,XNPTI,DISTI,XN1I,DN1I) 

C ICODE=26 

C WRITE ( IOUT, 1000 ) ICODE,XN II 

ZN10=XN1I(3) 

XN1I(3)=ZN10  + ZS(  1 ) 

HDIFF =HEIGHT(XN II, XRCI,XDB1,XDB2) 
IF(HDIFF.GT,20.0)GO  TO  113 

CALL  SECTN(XP1,XK,XRCI,XDB1,XDB2,XDB3,XDB4) 

C ICODE=27 
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on  no  ooo  no  on  oo 


C WRITE(IOUT,1000)ICODE,XDB3,XDB4 

CALL  NRPT(XDB3,XDB4,XRCl,XNPTJ,DISTJ ) 

CALL  NR1(XDB3,XDB4,XRCI,XNPTJ,DISTJ,XN2,DN2) 

C ICODE= 1 27 

C WRITE(IOUT,1000)ICODE,XNPTJ,XN2,DISTJ,DN2 

HDIFF=HEIGHT(XN2.XRCI,XRB1,XRB2) 

ICODE=227 

WRITE(IOUT,1002)ICODE,HDIFF 
IF(HDIFF,GT.-2.0)GO  TO  170 
DO  105  11=1, NQ 
IQ  = NQ+  1 -II 
XN1I(3)=ZN10+ZS(IQ) 
HDIFF=HEIGHT(XN1I,XRCI,XDB1,XDB2) 
DELP=DEL(XN1I,XRCI,XDB1»XDB2,HDIFF,DN1I) 
IC0DE=327 

WRITE (IOUT, 1002 ) ICODE,DELP 
IF(DELP.GE. 12.5,AND.IQ.EQ.2)GO  TO  170 
IF(DELP.GT,DELP1 (IQ) )G0  TO  100 
IF ( 10 . EQ . 1 )G0  TO  110 
GO  TO  105 
100  M D I F F = 1 

DELP1  ( IQ ) =DELP 
105  CONTINUE 
ICODE=28 

WRITE(IOUT,1000)ICODE,DELP1(1),DELP1(2) 

CALL  REPLCE ( XK,XR2  ) 

110  CONTINUE 
ICODE=29 

WRITE(IOUT,1000)ICODE,XR2 
DIFFRACTION  AFTER  REFLECTION 

CALL  IMAGE(XRB1,XRB2,XNPTI,XNPTJ) 

CALL  IMAGE(XRB1.XRB2,XR1,XR1I) 

115  CALL  IMAGE(XRB1 , XR82 , XR2 , XR2I ) 

ICODE=30 

WRITE (IOUT, 1000 )IC0DE,XR1X,XR2I,XNPTI 
IF(NRFDF,EQ,0)GO  TO  145 
DO  140  KRFDF=1 ,NRFDF 
KBAR2=KRDNUM(KRFDF) 

CALL  REPLCE(RB1 ( 1, 1 , KRFDF  ) ,XDB1  ) 

CALL  REPLCE(RB1 ( 1,2,KRFDF) ,XDB2) 

ICODE  = 3 1 

WRI TE (IOUT, 1001 )ICODE,KBAR2,XDB1,XDB2 
IF(KBAR2.E0.KBAR1 )G0  TO  140 
CALL  INTCPT(XRB1,XRB2,XRC,XR2I,XRB4) 

CALL  BLOCN(XRB3,XRB4,XRC,XDB1,XDB2,XJ,LOC) 

IF ( LOC . EQ , 0 ) GO  TO  140 
IF ( LOC  . EQ , 3 ) GO  TO  120 
CALL  INTCPT(XR1I,XR2I,XRC,XJ,XKI) 
XKI(3)=ZC0R(XR1I,XR2I,XKI) 

DELTA=-0.5 

CALL  MOVE(XKI,XKI,XR  1 1 , DELTA  , IERR  ) 

IF ( IERR  , EQ  , 4 ) RETURN 
I F ( LOC . N E , 1 ) GO  TO  135 
GO  TO  125 

120  CALL  REPLCE(XR2I,XKI) 

125  CALL  INTCPT(XR1I,XRC,XDB1,XDB2,XDB3) 
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C ICODE=32 

C WRITE ( IOUT, 1001 )ICODE,LOC,XJ,XKI,XDB3 

XDB3(3)=ZC0R(XDB1,XDB2,XDB3) 
HDIFF=HEIGHT(XR1I,XDB3,XRB1,XRB2) 

IF (HDIFF  .GT,-2.0)GO  TO  165 

CALL  NR1(XR1I,XKI,XRC,XNPTJ,DISTI,XN1I,DN1I) 
ZN  10  = XN  II ( 3 ) 

DO  130  11=1, NQ 

IQ=NQ+ 1 -II 

XN1I(3)=ZN10+ZS(IQ) 

HDIFF=HEIGHT(XN1I,XRC,XDB1,XDB2) 

IF (HDIFF. GT, 20.0. AND. IQ.EQ,2)GO  TO  140 
DELP=DEL(XN1I,XRC,XDB1,XDB2, HDIFF, DN1I) 

IF ( MDIFF.EQ, 1 ,AND,DELP.GT,-0,2 )GO  TO  165 
IF(DELP.LE.DELP2(I0) )GO  TO  140 
DELP2 ( IQ  ) =DELP 
130  CONTINUE 
C ICODE=33 

C WRITE(I0UT,  1000)ICODE,DELP2(  1 ) , DELP2 ( 2 ) 

IF(DELP2(2)  ,GE.  12.5)G0  TO  165 
135  CALL  REPLCE(XKI,XR2I) 

140  CONTINUE 

145  CALL  REPLCE(XR2I,XKI) 

IDXR=IDXR+1 
C ICODE=34 

C WRITE(IOUT,1001)ICODE,IDXR,XR2I 

IFfIDXR.LT. 11)QO  TO  150 
IERR=3 
RETURN 

150  DO  155  IQ  = 1 » NQ 

DELR(IQ,IDXR)=AMAX1(DELP1(IQ),DELP2(IQ)) 

155  CONTINUE 

DO  160  1=1,3 
X IMG ( I, IDXR ) =XRCI ( I ) 

160  CONTINUE 
GO  TO  165 

165  CALL  IMAGE ( XRB  1 , XRB2, XKI, XK ) 

170  CALL  REPLCE ( XK, XR2  ) 

C IC0DE=35 

C WRITE ( IOUT, 1 000 ) ICODE, XR2 

175  CONTINUE 
C ICODE=36 

C WRITE(I0UT, 1000)ICODE,XR1,XR2 

C BEGIN  BARRIER  FACTOR  COMPUTATION 
180  NIMG=IDXR 

ANG=ANGLE(XR1,XR2#XRC) 

ADST=ANG/DIST 

IF(KP0S(XNPT,XR2,XR1 ) ,EQ.  1 )G0  TO  190 
ANG2=ANG1-ANG 
GO  TO  195 
190  ANG2=ANG1+ANG 
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C CONTRIBUTION  FROM  DIRECT  R A I 
195  DO  205  IQ=1,NQ 

IF(NQS(MR,IQ) ,EQ.0)GO  TO  205 
DELP=DELP0 ( IQ ) 

DO  200  KF= 1 , NF 
FB(KF,IQ)=BARFAC(KF,DELP) 

200  CONTINUE 
ICODE=37 

WRITE (IOUT, 1002)ICODE,FB( 1 # IQ ) ,FB(9,XQ) 

205  CONTINUE 
C CONTRIBUTION  FROM  REFLECTIONS 
IF(NIMG.EQ,0)GO  TO  230 
DO  225  KIMG=1,NIMG 
DO  210  1=1,3 
XRCI(I)=XIMG(I,KIMG) 

210  CONTINUE 

ANGI=ANGLE(XR1 ,XR2,XRCI) 

CALL  NRPT(XR1,XR2#XRCI,XNPTI,DISTI) 

RATIO=(ANGI/DISTI)/ADST 

DO  220  IQ=1,NQ 

IF(NOS(MR,XQ) ,EQ.0)GO  TO  220 
DELP=DELR ( IQ, KIMG ) 

DO  215  KF=  1 , NF 

FB(KF,IQ)=FB(KF,IQ)+BARFAC(KF,DELP)*RATIO 

215  CONTINUE 
220  CONTINUE 
C ICODE=38 

C WRITE(IOUT,1002)ICODE*F3(1,1),FB(1,2)*FB(9,1),FB(9,2) 

225  CONTINUE 

C COMPUTE  MEAN  ENERGY  LEVEL 

230  CALL  NR1 (XR1 , XR2 , XRC , XN PT , DIST , X N 1 ,DN1 ) 

DO  232  IQ= 1 , NQ 
NQO=NQS(MR,IQ) 

IF(NOQ.EQ.0)GO  TO  232 
DO  232  1=1, NQQ 
XLA(I,IQ)=0 
232  CONTINUE 

DO  242  I F = 1 , N F 
IPbIF 

IF ( I F . EQ • 1 )IP  = 5 
A=10.**(-1,E-8*4,**IP*DN1) 

T 1=0. 

DO  240  IQ=  1 , NQ 

IF(NOS(MR,IQ) .EQ.0)GO  TO  240 
NQQ=NQS (MR, IQ) 

T 2 = 0 . 

DO  235  1=1, NQQ 

T3=XLREF(MR,I,IF,IQ)*A*FB(XF,IQ)*FG(IF) 

T2=T2+VEXPH (MR,I,IQ)*T3 
IF(NF.GT.6tAND.IF.EQ.  1 )GO  TO  235 
IF(NF.LE.6,AND.IF,NE. 1 )GO  TO  235 
XLA(I,IQ)=XLA(I,IQ)+T3 
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235  CONTINUE 

T1=T1+CQ(IQ)*T2 
240  CONTINUE 

XLE(IF)=XLE(IF)+ADST*T1 
242  CONTINUE 

SA  1 = SIN ( ANG1  ) 

S A 2 = S I N ( A N G 2 ) 

CA  1 =COS ( ANG  1 ) 

C A 2 = C 0 S ( ANG2 ) 

CPPEV=ANG 
DO  246  N = 1 # 4 
T 1 = 0. 

DO  245  IQ=1,NQ 

IF(NQS(MR,IQ) ,EQ.0)GO  TO  245 
NQQ=NQS(MR,IQ) 

T 2 = 0 . 

DO  244  1=1, NQQ 

T2=T2+VEXPH(MR,I,IQ)*XLA(I,IQ)**N 

244  CONTINUE 
T1=T1+CQ(IQ)**(N**2)*T2 

245  CONTINUE 
IT3  = 2*N-  1 
T3=IT3 

XKA(N)=XKA(N)+T1/DIST**IT3*CPREV 
A F N = N 

CPREV  = 0.5/AFN*(SA2*CA2**IT3~SA1*CA1**IT3  + T3«'CpREV) 

246  CONTINUE 

C ICODE=39 

C WRITE(IOUT,1002)ICODE,XLE(1),XLE(9),CAP2 

ANG  1 = ABS ( A NG  2 ) 

C ICODE=40 

C WRITE (I OUT, 1000)ICODE,XR2,XR2G,XR2D#XR20 

IF(DSQR(XR2,XR20).LT.1.)RETURN 
DELTA= 1 . 

CALL  MOVE(XR2,XR1,XR  10, DELTA, IERR) 

IF(IERR.EQ# 4)  RETURN 
C IC0DE=41 

C WRITE ( IOUT,  1000 ) ICODE,XR  1 

IF(DSQR(XR2,XB2G) ,LT.  1.0JGO  TO  247 
CALL  REPLCE ( XR2G, XR2  ) 

GO  TO  91 

247  IF ( DSQR ( XR2, XR2D ) ,LT, 1 ,0)GO  TO  25 
CALL  REPLCE(XR2D,XR2) 

GO  TO  72 

C1000  FORMAT ( 6H0CODE0I3, 6F9 , 2/6F9.2 ) 

Cl  001  FORMAT ( 6H0CODE®I3,I4, 6F9  #2/6F9,2  ) 

C 1 002  FORMAT (6H0CODBPI3#6E  12,2) 

END 

$ 1 
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C SUBROUTINES  FOR  TRAFFIC  NOISE  12/13/71 
SUBROUTINE  INTER(NR,IQ) 

DIMENSION  TL1  (2,9,2) ,TL(36) 

EQUIVALENCE (TL(  1 ) ,TL  ! ( 1#  1#  1 ) ) 

COMMON/BLK2/NQ 

COMMON/STORE4/XMPH ( 10,5, 2), VEXPH( 10,5,2) 

COMMON/DR  I V2/NQS ( 10,2) , N F 
COriMON/INTEB  1/XIREF  ( 10,5,9,2) 

DATA  TL ( 1 ) , TL(2),  TL(3),  TL(4),  TL(5),  TL ( 6 ) , TL(7),  TL(8), 

1 TL( 9 ) ,TL( 10) ,TL( 1 1 ) ,TL(  12) ,TL ( 13)  ,TL(  14  ) ,TL(  15)  ,TL(  16)  , 

2 TL(17),TL(18),TL(19),TL(20),TL(21),TL(22),TL(23),TL(24), 

3 TL(25) ,TL(26),TL(27) ,TL(28) ,TL(29)  ,TL(30)  ,TL(31  ),TL(32) , 

4 TL(33),TL(34),TL(35),TL(36)/ 

5 61. ,75. ,38. ,48. ,45. ,57. ,47. ,62, ,55. ,66, ,58 ,,70. ,54. ,72. ,49. 

6 42  . , 57  , , 

7 2*87., 2*60., 2*73. ,2*78. ,2*83,, 2*82, ,2*79, ,2*74. ,2*66,/ 

NQQ=NQS ( NR , IQ ) 

CONS=ALOG10(70./30. ) 

DO  10  I F = 1 , N F 

TEMP=TL1  (2,IF,IQ)-TL1 ( 1,IF,IQ) 

DO  10  1=1, NQQ 

XLPEF1=TL1 ( 1,IF,XQ)+TEMP*ALOG10(XMPH(NR,I,IQ)/30. )/CONS 
XLREF(NR,I,IF,IQ)  = 10.**( (XLREF1-66.  )/10.  ) 

10  CONTINUE 
RETURN 
END 
5 1 

FUNCTION  ACOS(X) 

IF(X.EQ»0,0)  GO  TO  5 
ACOS=ATAN(SQRT( 1 .0-X*X)/X) 

IF(X.LT.0.0)  ACOS=3. 1416+ACOS 
RETURN 

5 A C 0 S = 1 .5708 

RETURN 
END 

$1 

FUNCTION  AMAG ( X 1 V,X2V  ) 

C FIND  MAGNITUDE  OF  VECTOR 

DIMENSION  X1V(3),X2V(3) 

AMAG  = SQRT(DSQR(X1V,X2V)  ) 

RETURN 
END 
5 1 

FUNCTION  ANGLE(X1V,X2V,X3V) 

DIMENSION  X1V(3)  ,X2V(3)  ,X3V(3) 

D 13  = DSQR(X1V,X3V) 

D23=DSQR(X2V,X3V) 

ANGLE=  1 .5708 

IF (D13*D23,EQ,0. )RETURN 

D12=DSQR(X1V,X2V) 

ANGLE  = ACOS( ( D 2 3 + D 1 3 -D 1 2 ) / ( SQRT ( D 1 3 * D 23  ) * 2 , ) ) 

RETURN 

END 


,63,, 
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FUNCTION  BAP.FAC(KF,DELP) 

C FIND  BARRIER  FACTOR 

IF (DELP,EQ,-0.2)GO  TO  3 
IF (DELP.GE, 1 2 , 5 ) GO  TO  4 
IP  = KF 

IF(KF.EQ. 1 )IP=5 

A=DELP*(2.**IP)/5.7 

IF(A.GE.78.5)GO  TO  4 

IF( A.GT.0. ,AND.A.LT.78.5)GO  TO  5 

IF( A.EO.0. ) GO  TO  6 

IF( A.GT.-1 .25. AND. A.LT.0. )GO  TO  7 

3 BARF  AC=  1.0 
RETURN 

4 BARFAC=4 , 0E-3 
RETURN 

5 BARFAC=(TANH(SQRT(A) )**2)/A/3. 16 
RETURN 

6 BARFAC=.316 
RETURN 

7 A 1 = AB S ( A ) 

BARFAC=(TAN(SQRT(A1  ) ) )**2/A1/3t  16 

RETURN 

END 

% 1 

SUBROUTINE  BLOCN(X1V,X2V,X3V,X4V,X5V,X6V,LOC) 

C FIND  RELATIVE  LOCATION  OF  BARRIER 

DIMENSION  X1V(2),X2V(2),X3V(2)eX4V(2),X5V(2),X6V(2),XAV(2),XBV(2) 
LOC=KCUT(X1V,X3V,X4V,X5V) 

LOC=LOC+KCUT(X2V,X3V,X4 V,X5V) *2 
IF(LOC.E0.3)RETURN 

CALL  TRI(X1V,X2V,X3V,X4V,XAV,KTRI) 

IF ( LOC . EQ , 0 ) GO  TO  4 
IF ( KTRI . EQ , 1 ) GO  TO  6 

2 CALL  INTCPT(X1V,X2V#X3V,X5V,XBV) 

IF ( LOC  , EQ  . 4 ) GO  TO  5 

3 X 6 V ( 1 ) =XB V ( 1 ) 

X6V(2)=XBV(2) 

RETURN 

4 IF(KTRI.EQ,0)RETURN 
LOC  = 4 

GO  TO  2 

5 IF(KPOS(X1V,XAV,XBV) .EQ. 1 )GO  TO  3 

6 X 6 V ( 1 ) =X A V ( 1 ) 

X 6 V ( 2 ) =X  A V ( 2 ) 

RETURN 

END 

$ 1 

FUNCTION  DEL(X1V,X2V,X3V,X4V,HDIFF,DN1) 
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c find  path  length  difference 

DIMENSION  X1V(3).X2V(3),X3V(3),X4V(3),XAV(3),XBV(3) 

CALL  NRPT(X3V,X4V,XlV,XAV,DISTA) 

CALL  NRPT(X3V,X4V,X2V,XBV,DISTB) 

DISTCaDSQR (XBV,XAV) 

DEL=SQRT ((DISTA+DISTB)**2+DISTC)-DN1 

IF(HDIFF.GT.0,)DEL=~DEL 

RETURN 

END 

$ 1 

FUNCTION  DSQR ( XI V,X2V ) 

DIMENSION  X1V(3),X2V(3) 

DSQR=0.0 
DO  10  1=1,3 

10  DSQR=(X1V(I)-X2V(I)  )**24'DSQR 
RETURN 
END 
S 1 

SUBROUTINE  ENDPT ( X 1 V , X2 V , X 3 V . X4 V , X5 V , X6 V , KTRIG , IERR ) 

c find  new  endpoint 

DIMENSION  X 1 V ( 3 ) , X2V ( 3 ) , X3 V ( 3 ) , X4V ( 3 ) , X5V ( 3 ) , X6V ( 3 ) , XDUM ( 3 ) 

IERR=0 

KTRIG=0 

ITRIG= 1 

CALL  REPLCE ( X 1 V, XDUM ) 

1 CALL  BLOCN(XDUM,X2V,X3V,X4V,X5V,X6V,LOC) 

IF ( LOC . EQ , 0 ) RETURN 

IF(LOC.NE,3)  GO  TO  2 
CALL  REPLCS(X2V,X6V) 

GO  TO  4 

2 X6V(3)=ZCOR(X1V,X2V,X6V) 

IF ( ITRIG  • EQ , 0 ) GO  TO  5 
IF(AMAG(X1V,X6V) .GT.0.51 ) GO  TO  5 
I T R I G = 0 

DELTA=-0,51 

CALL  MOVE (XDUM, XDUM, X2V, DELTA, IERR) 

GO  TO  1 

5*  DELTA=-0,5 

IF ( LOC  , EQ  . 1 ) GO  TO  3 

CALL  MOVE (X6V,X2V, XIV, DELTA, IERR) 

RETURN 

3 CALL  MOVE(X6V,X6V, XIV, DELTA, IERR) 

4 KTR IG= 1 
RETURN 
END 

SI 

FUNCTION  GAMMA(X) 

DATA  Cl/2,5066272/ 

H = X 
P=1  . 

IF ( H , GT , 25 , ) G 0 TO  15 
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5 IF(H.GE,6,  )G0  TO  10 
P = P/H 
H = H+  1 , 

GO  TO  5 
10  Y=1./H 

GAMMA=EXP  ( ..  H ) * H * * ( H - 0 , 5 ) * C 1 * P 
12  GANMA*GAMMA* ( 1 .♦( 1 .+Y/24. ) *Y/12, ) 

RETURN 

15  TEMP=-H+(H-0.5) *ALOG(H) 

GAMMA  = C 1 *EXP ( TEMP ) 

GO  TO  12 
RETURN 
END 
S 1 

FUNCTION  GAUSS ( Y ) 

C GAUSSIAN  PROBABILITY  FUNCTION 
X = ABS ( Y ) 

G A U S S = 0 . 

IFf X.GT.5. )GO  TO  10 
T=1 ./( 1 .+0.33267*X) 

GAUSS  = 0 ,436  1836 *T-0*  1 20  1676 *T*T  + 0 , 93729b *T*T*T 
GAUSS =GAUSS/EXP(X*X/2.)/2, 50663 
10  IF( Y.GT.0, )GAUSS=1 , -GAUSS 
RETURN 
END 

* 1 

FUNCTION  HEIGHT (X IV, X2V,X3V,X4V) 

C FIND  HEIGHT  DIFFERENCE 

DIMENSION  X1V(3),X2V(3),X3V{3),X4V(3),XI(3) 

CALL  INTCPT(X1V,X2V,X3V,X4V,XI) 

HEIGH T=ZCOR (XIV, X2V,XI)-ZCOR(X3V,X4V, XI) 

RETURN 
END 
5 1 

FUNCTION  IEPS  (XIV,  X2V,  X3V,  X4V,  DELI) 

C CHECK  ON  PATH  LENGTH  DIFFERENCE 

DIMENSION  X1V(3),X2V(3)#X3V(3),X4V(3) 

IEPS  =0 

DIST  = AMAG  (XIV, X2V) 

HDIFF  = HEIGHT  ( X 1 V , X 2 V , X 3 V , X 4 V ) 

DEL  2 = DEL  (XIV,  X2 V , X3V,  X4V,  HDIFF,  DIST) 

DELM=  (DELI +DEL2 ) /2  • 

IF(  (ABS(DEL2-DEL1)  - 0.1 -DELM/50 . * (1.+DELM))  ,GT,0,)  IEPS 

RETURN 
END 
SI 

SUBROUTINE  IMAGE ( XIV, X2V,X3V,X4V) 
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C FIND  IMAGE  POINT 

DIMENSION  X1V(3),X2V(3),X3V(3),X4V(3) 

AX  = X2V ( 1 ) -X  1 V ( 1 ) 

AY=X2V(2)-X1V(2) 

AXY=AX**2+AY**2 
R ATIO=0 . 

IF ( AXY  . EQ  « 0 , ) GO  TO  10 

RATIO=(((X3V(2)"X2V(2))*AX~(X3V(1)-X2V(1))*AY)*2.0)/AXY 
10  X4V ( 1 )=X3V(  1 )+AY*RATIO 
X4V(2)=X3V(2)-AX*RAHO 
X4V(3)=X3V(3) 

RETURN 
END 
S 1 

SUBROUTINE  INTCPT(XlV,X2V#X3V,X4V,X5V) 

C FIND  INTERCEPT  OF  TWO  LINES  IN  A PLANE 

DIMENSION  X1V(2),X2V(2),X3V(2).X4V(2),X5V(2) 

A X = X 2 V ( 1 )-X1V(  1 ) 

AY=X2V ( 2 ) -X  1 V ( 2 ) 

BX  = X4 V ( 1 ) -X3V  ( 1 ) 

B Y = X 4 V ( 2 ) -X3V ( 2 ) 

C 1=AY*X2V( 1 )-AX*X2V(2) 

C2=BY*X4V ( 1 ) -BX*X4V ( 2 ) 

D=AX*BY-AY*BX 
IF(D**2,LT,1.E-6)GO  TO  2 
X 5 V ( 1 ) = (AX*C2-BX*C1 )/D 
X 5 V (2)  = (AY*C2-BY*C1 )/D 
RETURN 

2 D=SQRT(AX**2+AY**2) 

X5V ( 1 )=X1V { 1 ) + ( AX/D ) * 1 .E+15 

X5V(2)=X1V(2)+(AY/D)*1.E+15 

RETURN 

END 

$ 1 

FUNCTION  KCUT(X1V,X2V,X3V,X4V) 

C DETERMINE  IF  TWO  LINE  SEGMENTS  CROSS 

DIMENSION  X1V(2).X2V(2),X3V(2)#X4V(2),X5V(2) 

K C U T = 0 

CALL  INTCPT(X1V,X2V,X3V,X4V,X5V) 

IF(KPOS(X1V,X2V#X5V).NE.1)RETURN 
IF(KPOS(X3V,X4V,X5V).EQ,1)KCUT=1 
RETURN 
END 
S 1 

FUNCTION  KP0S(X1V,X2V,X3V) 

C FIND  POSITION  OF  POINT  ON  LINE 

DIMENSION  X1V(2) ,X2V(2) »X3V(2) 

K PO  S = 1 

IF(((X3V(1)«X1V(1))*(X3V(1)«X2V(1))+(X3V(2)SX1V(2))*(X3V(2)-X2V(2) 
1 ) ) . GT , 0 . )KPOS=0 
RETURN 
END 
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S 1 

SUBROUTINE  MIDP  (XIV,  X2V,  X3V ) 

C FIND  CENTER  POINT 

DIMENSION  X 1 V ( 3 ) , X2V(3),  X3V(3) 

DO  10  1=1,3 

10  X 3 V ( I ) = (X1V(I)+X2V(I)  )/2, 

RETURN 

END 

% 1 

SUBROUTINE  M0VE(X1V,X2V,X3V,DELTA,IERR) 

C MOVE  ENDPOINT  OF  ROAD 

DIMENSION  X1V(3),X2V(3),X3V(3) 

IERR  = 0 

TEMP=AMAG ( X 1 V, X3V ) 

IF ( TEMP , EQ , 0 . ) GO  TO  3 

FCTR=DELTA/TEMP 

DO  2 1=1,3 

X2V(I)=X1V(I)+(X1V(I)-X3V(I) )*FCTR 

2 CONTINUE 
RETURN 

3 IERR=4 
RETURN 
END 

$ 1 

SUBROUTINE  NPPT(X1V,X2V,X3V,X4V,DIST) 

C FIND  NEAREST  POINT  TO  LINE 

DIMENSION  X1V(3),X2V(3),X3V(3),X4V(3),AV(3),BV(3) 

EQUIVALENCE  (AV(1),AX),(AV(2),AY),(AV(3),AZ),(BV(1),BX),(BV(2),BY 
1 ) , (BV(3) ,BZ) 

DO  5 1=1,3 
AV (I)=X2V ( I ) -X 1 V (I) 

BV(I)=X3V(I)-X1V(I) 

5 CONTINUE 
RATIO=0. 

TEMP  = DSQR ( X 2 V , X IV ) 

IF(TEMP.NE,0. ) RATIO* ( AX*BX+A Y*BY+AZ*BZ ) /TEMP 
DO  10  1=1,3 

X4V ( I ) =X 1 V (I)+RATIO*AV (I) 

10  CONTINUE 

DIST=AMAG(X4V,X3V) 

IF(DIST.EQ,0.  )DIST=1  . 

RETURN 
END 
5 1 

SUBROUTINE  N R 1 ( X 1 V , X 2 V , X 3 V , X 4 V , D I S T , X 5 V , D N 1 ) 
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C FIND  NEAREST  POINT  TO  LINE  SEGMENT 

DIMENSION  X1V(3),X2V(3)#X3V(3),X4V(3),X5V(3) 
IF(KPOS(X4V,X2V,X1V) .EQ, 1 )GO  TO  2 
IF(KPOS(X1V,X4V,X2V) ,EQ. 1 )GO  TO  4 
CALL  REPLCE(X4V,X5V) 

DN  1 =DIST 
RETURN 

2 CALL  REPLCE(X1V,X5V) 

GO  TO  6 

4 CALL  REPLCE(X2V,X5V) 

6 DN1=AMAG(X5V,X3V) 

RETURN 
END 
S 1 

SUBROUTINE  REPLCE ( X 1 V , X2V ) 

DIMENSION  X1V(  3)  ,X2V(3) 

X 2 V ( 1 ) =X 1 V ( 1 ) 

X 2 V ( 2 ) =X 1 V ( 2 ) 

X2V(3)=X1V(3) 

RETURN 

END 

% 1 

SUBROUTINE  SECTN(X1V,X2V,X3V,X4V,X5V,X6V,X7V) 

C FIND  EFFECTIVE  BARRIER  SECTION 

DIMENSION  X1V(3),X2V(3),X3V(3),X4V(3),XbV(3),X6V(3),X'7V(3) 
CALL  INTCPT(X1V,X3V,X4V,X5V,X6V) 

XbV(3)=ZC0R(X4V,X5V,X6V) 

CALL  INTCPT(X2V,X3V,X4V,X5V,X7V) 

X7V(3)=ZCOR(X4V,X5V,X7V) 

RETURN 
END 
S 'I 

FUNCTION  TAN ( X ) 

TAN=SIN(X)/COS(X) 

RETURN 

END 

$ 1 

SUBROUTINE  TRI(X1V,X2V,X3V,X4V,X5V,KTRI) 

C FIND  IF  POINT  IN  TRIANGLE  AND  LOCATE  INTERCEPT 

DIMENSION  X1V(2),X2V(2)»X3V(2),X4V(2)#X5V(2) 

CALL  INTCPT(X1V,X2V,X3V,X4V,X5V) 

KTRI=0 

IF(KPOS (XIV, X2V,X5V).EQ,0)  RETURN 

IF(KPOS(X3V,X5V,X4V).EQ.1)KTFI=1 

RETURN 

END 

$ 1 

FUNCTION  XINGAM ( X # TM ) 
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C INCOMPLETE  GAMMA  FUNCTION 
TTM=TM 

IF  (TTM.LT.50)  GO  TO  10 

XINGAM=1,-GAUSS(SQRT(2.*X)..SQRT(4,*TTM-1,  ) ) 

RETURN 
10  XINGAM=0 

IF  (TTM.GT.15.)  GO  TO  30 

TEMP=(X/2,  )**TTM/GAMMA(TTM+1 ,)/EXP(X/2.) 

20  XINGAM=XINGAM-TEMP 

IF (TEMP.LT, 1 . E-20 ) GO  TO  40 
TTM=TTM+1 . 

IF ( TTM.GT.15. ) GO  TO  30 
TEMP=TEMP*X/2./TTM 
GO  TO  20 

3 0 T E M P = ((X/2,/TTM)*+(1./3.)-1,  + 1t/9,/TTM)*SQRT(9,*TTM) 
XINGAMsXlNGAM+l  .-.GAUSS(TEMP) 

RETURN 

40  XINGAM=XINGAM+1 , 

RETURN 

END 

FUNCTION  XPROB ( A2, AA3, AA4, A5, X ) 

C CALCULATION  OF  PROBABILITY 

P=1 . -XINGAM ( 2*X, AA3 ) 

F 1 = A 5 * ( 1 . -A2 ) *AA4 
F2=X/AA3 
N = 0 
S 2=  1 . 

10  Q=P-F1*S2 

IF(F1*S2/Q.LE. 0.001 )GO  TO  20 
N = N+  1 
X N = N 

F2=F2*X/( AA3+XN ) 

S2=S2+F2 

F1«F1*(1.-A2)*(AA4  + XN)/(XN+1,  ) 

P = Q 

GO  TO  10 
20  X P R 0 B = Q 
RETURN 
END 


FUNCTION  ZCOR(X1V,X2V,X3V) 

C FIND  Z COORDINATE 

DIMENSION  X 1 V ( 3)  ,X2V(3) #X3V(3) 

TEM  1 =X2 V ( 1 )-X1V( 1 ) 

TEM2=X2V(2)»X1V(2) 

TEM3  = X2V(3)t*X1V(3) 

IF(ABS(TEM1  ) .GT,  ABS ( TEM2 ) )GO  TO  10 
ZCOR=X1V ( 3)+(X3V ( 2 ) -X 1 V ( 2) ) ♦TEM3/TEM2 
RETURN 

10  ZCOR  = X1V(3)  + (X3V(1)-X  i V ( 1 ) ) *TEM3/TEM 1 
RETURN 
END 


5 1 
S0 
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4.  USER'S  GUIDE  TO  THE  COMPUTER  PROGRAM 


4.1  Format  and  Sequence  of  Data  Input 

The  computer  program  is  written  to  accept  input  data  from  a card 
reader.  Three  types  of  input  data  are  allowed: 

Integer  — a fixed-point  number  written  without  a decimal  point. 
All  integers  must  be  right-adjusted  within  the  allotted  field. 

Real  Constant  — a floating  point  number  written  with  a decimal 
point.  The  real  constant  may  be  situated  anywhere  within  its 
allotted  columns. 

Alphanumeric  — any  combination  of  alphabetic  and  numeric  char- 
acters. Alphanumeric  data  may  be  situated  anywhere  within  its 
allotted  columns. 

The  first  card  read  by  the  program  is  a title  card  with  arbi- 
trary alphanumeric  information  that  may  appear  in  columns  2 
through  60.  Up  to  five  blocks  of  data  may  follow,  in  arbitrary 
sequence,  with  information  about  the  following: 

1 Program  initialization  parameters, 

2 Road  and  vehicle  parameters,  . 

3 Barrier  parameters, 

4 Ground  cover  parameters, 

5 Receiver  data. 

Each  data  block  must  be  preceded  by  a control  card  containing 
an  integer  between  1 and  5 in  column  5-  This  integer  char- 
acterizes the  data  block  as  listed  above.  A card  with  the 
integer  6 in  column  5 indicates  the  end  of  the  input  data 
string . 
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Multiple  cases  may  be  run  back  to  back,  whereby  each  case  starts 
with  a title  card  and  terminates  with  a card  having  a 6 in 
column  5.  The  entire  situation  need  not  be  redefined  for  each 
additional  case;  rather,  it  is  sufficient  to  specify  only  those 
aspects  of  the  problem  that  differ  from  the  case  immediately 
preceding. 

The  first  of  a series  of  cases  must  specify  initialization  para- 
meters, vehicle  and  road  parameters,  and  receiver  locations. 

If  no  information  is  provided  for  barriers  for  the  first  case, 
the  computer  program  assumes  that  there  are  no  barriers,  and 
all  calculations  related  to  barrier  diffraction  and  reflection 
are  bypassed.  Ground  absorption  is  treated  in  the  same  manner. 
Once  barrier  information  has  been  provided,  however,  failure 
to  specify  barrier  data  in  subsequent  cases  of  the  same  series 
of  cases  will  not  reinitialize  the  number  of  barriers  to  zero. 
Instead,  the  barrier  data  appropriate  to  the  preceding  case 
will  be  retained. 

Specifications  for  the  five  types  of  data  blocks  accepted  by 
the  program  are  given  below. 


4.1.1  Program  initialization  parameters 

Initialization  parameters  must  appear  in  the  first  of  a series 
of  cases.  A control  card  is  provided  with  the  integer  1 in 
column  5 to  identify  the  data  block.  Six  data  cards  must 
follow,  in  arbitrary  sequence,  to  provide  information  about: 

1 Receiver  height  adjustment, 

2 Number  of  frequency  bands , 

3 Standard  deviation  for  cars. 
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4 Standard  deviation  for  trucks, 

5 Height  adjustment  for  cars, 

6 Height  adjustment  for  trucks. 

An  initialization  parameter  is  entered  as  a real  constant  in 
columns  1 through  10  of  each  data  card.  Column  15  contains  an 
integer  between  1 and  6 to  identify  the  parameter,  as  indicated 
in  the  above  list,  and  alphanumeric  information  is  contained 
in  columns  31  through  60.  The  last  of  the  set  of  six  cards 
must  contain  the  letter  L in  column  20  to  signify  the  end  of 
the  data  block. 

The  receiver  height  adjustment  is  a height  in  feet  to  be  added 
to  the  z-coordinate  of  all  receivers  specified  in  data  block  5. 
For  normal  conditions,  it  will  be  0.0.  For  special  cases, 
such  as  the  comparison  of  levels  at  the  ground  floor  of  a 
building  with  levels  at  a higher  floor,  a height  adjustment 
of  10.0  feet  per  floor  is  recommended.  Negative  height  adjust- 
ments are  allowable.  Such  an  adjustment  will  not  cause  the 
program  to  assume  that  the  receiver  is  shielded  by  the  ground. 

The  number  of  frequency  bands  for  which  calculations  are  to  be 
performed  is  indicated  by  a number  (entered  as  a real  constant) 
between  1.0  and  9*0.  By  using  the  number  1.0,  all  attenuations 
due  to  atmospheric  absorption,  ground  absorption,  and  barrier 
diffraction  will  be  calculated  for  a frequency  of  500  Hz. 

This  frequency  is  most  typical  for  the  spectrum  of  road  traffic 
noise.  By  using  the  number  9*0,  the  attenuations  will  be  cal- 
culated in  addition  in  eight  octave  bands  with  center  frequen- 
cies ranging  from  62.5  Hz  to  8000  Hz.  The  final  results  will 
be  more  accurate  in  this  case,  but  total  computation  time  will 
be  increased.  By  using  any  number  n between  1.0  and  9-0,  the 
attenuations  will  be  calculated  for  500  Hz  and  for  n-1  octave 
band  starting  from  the  62.5  Hz  band. 
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The  standard  deviation  for  cars  and  trucks  is  the  standard 

deviation  of  reference  levels  at  a short  distance  (r  = 50  ft) 

o 

from  individual  vehicles.  The  value  recommended  for  cars  is 
2.5  (dB)  and  for  trucks  3-5  (dB). 

The  height  adjustment  for  cars  should  be  0.0,  since  tire  noise 
is  assumed  to  be  the  dominant  noise  source.  Other  height  ad- 
justments for  engine  or  fan  noise  should  not  be  made. 

The  height  adjustment  for  trucks  is  typically  8.0  (ft)  for  ex- 
haust noise.  Other  height  adjustments  for  engine  or  fan  noise 
will  result  in  slightly  higher  attenuation  by  barrier  diffrac- 
tion. A height  adjustment  0.0  for  tire  noise  is  incompatible 
with  the  noise  spectrum  stored  in  the  program. 

4.1.2  Road  and  vehicle  data 

The  control  card  for  this  data  block  must  have  the  integer  2 in 
column  5 and  a second  integer  right-adjusted  in  columns  6 through 
10  to  indicate  the  number  of  roadways.  A maximum  number  of  10 
roadways  is  allowed. 

The  number  of  roads  to  be  specified  depends  on  the  number  of 
traffic  arteries  with  constant  hourly  traffic  flow  that  one 
wants  to  consider.  For  a receiver  located  far  from  a multi- 
lane highway  without  ramps,  consideration  of  a single  traffic 
artery  is  sufficient.  For  receiver  locations  closer  to  the 
highway,  many  traffic  arteries  in  parallel  might  be  considered 
to  account  for  individual  lanes.  Ramps  to  highways  and  inde- 
pendent highways  are  treated  as  separate  roads. 

Following  the  control  card,  a set  of  one  or  more  cards  is 
required  to  provide  vehicle  data  for  the  first  roadway . 
card  contains  a real  constant  in  columns  1 through  10  to  indi- 
cate vehicles  per  hour,  a real  constant  in  columns  11  through 
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20  to  indicate  vehicle  speed  in  miles  per  hour,  and  the  integer 
1 or  2 in  column  25  to  indicate  vehicle  type.  The  final  card 
of  this  set  has,  in  addition  to  the  above  information,  an  L 
in  column  31. 

An  integer  1 in  column  25  stands  for  passenger  cars,  a 2 for 
trucks.  Up  to  five  cards  for  each  vehicle  type  can  be  used  in 
random  order  to  specify  hourly  traffic  flow  at  different  average 
speeds.  Unless  the  speed  distribution  is  known  from  measure- 
ments, the  "Highway  Capacity  Manual"  is  recommended  for  deter- 
mining the  average  speed  from  the  hourly  traffic  volume  and 
the  average  highway  speed. 

Following  the  vehicle  information  for  the  first  road  are  cards 
with  Cartesian  coordinates  for  the  end  points  of  straight-line 
approximations  to  the  first  road.  The  x-y  plane  is  parallel  to 
sea  level  and  the  z-coordinate  gives  the  height  above  sea  level 
or  above  any  plane  parallel  to  sea  level.  All  coordinates  are 
in  feet.  The  x,  y,  and  z-coordinates  are  entered  as  real 
constants  in  the  fields  bounded  by  columns  1 and  10,  11  and 
20,  and  21  and  30,  respectively.  The  card  with  the  last  point 
of  the  road  has  an  L in  column  31.  A roadway  may  contain  up 
to  eleven  end  points  (i.e.,  ten  sections).  A road  containing 
more  than  ten  sections  should  be  treated  as  two  or  more  roads, 
each  having  fewer  than  ten  sections. 

Vehicle  and  road  data  are  specified  for  the  remaining  roadways 
in  a similar  manner.  Data  must  be  specified  for  as  many  roads 
as  indicated  by  the  integer  contained  in  the  second  field  of 
the  control  card. 
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4.1.3  Barrier  parameters 

The  data  block  with  Information  about  obstacles  (i.e.,  "barriers") 
in  the  sound  propagation  paths  is  headed  by  a card  having  a 3 in 
column  5 and  an  integer  right-adjusted  in  columns  6 through  10 
to  indicate  the  number  of  barriers.  A maximum  number  of  10 
barriers  is  allowed. 

The  top  contour  of  a barrier  is  approximated  by  a straight-line 
segment,  and  no  sound  is  assumed  to  penetrate  below  this  con- 
tour. A single  barrier  may  contain  up  to  ten  sections.  The 
end  points  of  these  sections  are  specified  in  the  same  format 
as  the  end  points  of  road  sections,  except  that  the  last  point 
of  a line  is  identified  by  an  A or  an  R in  column  31  (rather 
than  L as  for  roads). 

An  R in  column  31  indicates  that  the  preceding  points  describe 
the  top  line  of  a rigid  plane  oriented  perpendicularly  to  the 
ground,  such  as  artificial  barriers  without  absorbing  material, 
facades  of  buildings,  rigid  walls  of  a depressed  highway,  etc. 

An  A in  column  31  indicates  that  the  preceding  points  describe 
the  top  line  of  a tilted  barrier,  a barrier  with  absorbing 
material,  an  earth  berm,  a hill,  or  some  other  obstacle  that 
reflects  sound  either  weakly  or  towards  the  sky,  directly  or 
via  a ground  reflection. 

The  top  line  of  a barrier  must  not  cross  a road. 

4.1.4  Ground  cover  parameters 

The  control  card  for  this  data  block  contains  a 4 in  column  ? 
and  an  integer  in  columns  6 through  10  to  indicate  the  number 
of  absorptive  ground  strips.  A maximum  number  of  10  ground  strips 
is  allowed. 
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The  areas  of  ground  cover  are  described  by  the  center  line  and 
the  width  of  rectangles.  The  x,  y,  and  z- coordinates  of  one 
end  point  of  the  center  line  are  given  as  real  constants  in 
the  fields  between  columns  1 and  10,  11  and  20,  and  21  and  30, 
respectively  of  a single  card.  The  same  card  contains  the 
width  of  the  rectangle  in  real  formats  in  the  field  between 
columns  31  and  40.  The  x,  y,  and  z-coordinates  of  the  other 
end  point  of  the  center  line  are  written  on  the  next  card 
together  with  a G or  a T in  column  31* 

A G identifies  the  ground  cover  as  high  grass  or  shrubbery, 
and  a T identifies  the  ground  cover  as  trees. 

The  rectangles  with  ground  cover  must  not  cross  a road. 

4.1.5  Receiver  data 

The  control  card  contains  a 5 in  column  5 and  an  integer  in 
columns  6 through  10  to  indicate  the  number  of  receivers.  A 
card  is  provided  for  each  receiver  to  indicate  the  x,  y,  and 
z-coordinates  of  the  receiver  location.  As  an  option,  background 
can  be  specified  for  each  receiver  in  terms  of  the  mean  and  the 
standard  deviation  of  normally  distributed  A-weighted  sound  pres- 
sure levels.  The  coordinates  are  entered  as  real  constants  in 
the  fields  bounded  by  columns  1 and  10,  11  and  20,  and  21  and  30, 
respectively.  The  mean  and  standard  deviation  of  the  background 
noise  level  can  be  entered  in  the  same  card  in  the  fields  bounded 
by  columns  31  and  40,  and  4l  and  50,  respectively. 

A maximum  of  fifteen  receivers  is  allowed. 

A receiver  cannot  be  located  on  a road,  nor  on,  over  or  under- 
neath a top  line  of  an  obstacle  (barrier),  nor  on  a ground  strip. 
However,  a receiver  can  be  located  between  two  adjacent  ground 
strips  if  the  location  is  not  identical  with  an  end  point  of  the 
ground  strip  center  line. 
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in  the  following  cases. 


4 . 2 Error 

Errors  are 

a.  If  the 
comput 


Messages 

detected  by 

top  line  of 
r prints  out 


the  computer 
an  obstacle  intersects 


a roadway. 


the 


ILLEGAL  BARRIER  INTERSECTS  ROADWAY 


together  with  the  following  data:  receiver  number,  roadway 

number,  road  section  number,  barrier  number,  barrier  section 
number . 


The  computer  then  proceeds  with  the  next  case. 

b.  If  the  center  line  of  a ground  strip  intersects  a roadway , 
the  computer  prints  out 

ILLEGAL  GROUND  STRIP  INTERSECTS  ROADWAY 

together  with  the  following  data:  receiver  number,  roadway 

number,  road  section  number,  ground  strip  number. 

The  computer  then  proceeds  with  the  next  case. 

c.  If  the  number  of  reflections  contributing  to  the  sound  level 
at  a receiver  exceeds  10,  the  computer  prints  out 

TOO  MANY  REFLECTIONS 

together  with  the  following  data:  receiver  number,  roadway 

number,  road  section  number. 

The  computer  then  proceeds  with  the  next  receiver. 

d.  Should  the  geometry  of  the  situation  be  such  that  a line 
segment  that  is  needed  for  computation  appears  to  have 
length,  the  program  continues  with  the  next  receiver  wirn- 
out  giving  an  error  message. 
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4.3  Data  Output 


The  data  output  starts  immediately  after  the  input  file  is  read. 
The  heading  TRAFFIC  NOISE  PREDICTION  is  typed  and  input  data 
are  then  typed  in  the  following  order: 

Title  card; 

Program  initialization  parameters  (for  the  first  of  a series  of 
cases  only ) ; 

Vehicle  and  road  parameters,  where  all  parameters  of  Type  1 
vehicles  (passenger  cars)  on  a given  road  are  typed  first, 
following  by  parameters  for  Type  2 vehicles  (trucks); 

Barrier  parameters,  with  the  type  of  barrier,  A or  R,  given 
in  parentheses  after  the  barrier  number; 

Ground  cover  parameters,  with  identification  G or  T following 
the  ground  strip  number  in  parentheses; 

Receiver  data. 

The  title  card  is  then  printed  out  again  together  with  a head- 
line for  the  receiver  identification.  The  output  consists  of 
the  following  data  for  each  receiver.  After  the  energy  mean  level 
for  a given  receiver  has  been  calculated,  the  receiver  number  and 
coordinates  are  typed.  If  a calculation  in  octave  bands  is  per- 
formed, the  octave  band  center  frequencies  are  typed  out;  under- 
neath them  are  the  calculated  A-weighted  octave  band  levels.  The 
final  typeouts  consist  of  the  A-weighted  overall  level  of  the 
energy  mean,  LE(A),  the  noise  pollution  level,  LNP,  and  percentile 
levels,  L90,  L50,  and  L10.  If  a background  noise  is  specified 
for  the  receiver,  the  last  line  gives  the  increase  in  LE(A),  in 
LNP,  and  in  the  percentile  levels  above  the  background  noise. 
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4.4  Recommendations  for  Usage 

The  computation  time  required  to  compute  sound  levels  at  a given 
receiver  increases  with  the  number  of  roadways,  barriers,  and 
ground  strips,  but  it  is  independent  of  the  total  number  of 
receivers.  Therefore,  if  one  wishes  to  analyze  a situation 
in  which  many  receivers  and  barriers  are  located  along  or 
near  a highway  many  miles  in  length,  we  suggest  that  the  prob- 
lem be  modularized  so  that  a smaller  geographical  area  can  be 
associated  with  each  receiver. 

Maps  showing  the  elevation  of  the  terrain  should  be  used  for  a 
detailed  description  of  the  roadway  and  of  top  lines  of  hills 
(barriers  of  Type  A) . It  should  be  borne  in  mind  that  the 
computer  program  accounts  for  5 dB  attenuation  for  sound  rays 
grazing  over  hills. 

Results  obtained  when  ground  absorption  is  considered  should  be 
interpreted  with  caution.  Because  of  the  lack  of  reliable  data 
in  this  area,  results  from  the  first-order  approximations  used 
in  the  computer  program  are  intended  for  comparison  with  field 
data  rather  than  for  prediction  purposes. 

4 . 5 Test  Cases 

As  an  example,  a listing  of  input  data  for  two  cases  and  the 
corresponding  computer  printout  is  given  on  the  following 
pages.  A graph  of  the  geometry  in  the  x-y  plane  is  shown  on 
Fig.  13. 

In  the  example,  test  case  1 involves  two  roadways,  a -ree  -o:.e, 
and  two  receivers.  The  curved,  major  roadway  with  10CR  passen- 
ger cars  per  hour  at  an  average  speed  of  6 0 mph  and  with.  -0^. 
trucks  per  hour  at  50  mph  is  approximated  by  two  st  . 
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FIG. 13  PLAN  VIEW  OF  THE  GEOMETRY  USED  IN  THE  TEST  CASES 
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sections.  The  straight,  minor  roadway  with  300  passenger 
cars  per  hour  at  an  average  speed  of  40  mph  has  only  one 
section . 

Since  most  of  the  sound  is  radiated  from  the  major  roadway  and 
the  receiver  locations  are  symmetrical  about  its  first  section, 
the  sound  pressure  levels  calculated  at  the  two  receivers  are 
about  equal. 

Test  case  2 involves  the  same  parameters  as  in  test  case  1 as 
well  as  two  barriers.  One  of  the  barriers  might  be  a hill  or 
an  earth  berm  and  is  therefore  of  an  absorptive  type.  The 
other  barrier  might  be  a building  and  is  specified  by  four 
reflective  sections. 

The  absorptive  barrier  is  between  10  and  20  ft  higher  than  the 
major  roadway  and  shields  receiver  no.  1 against  noise  from 
almost  the  entire  major  roadway.  However,  diffraction  of 
sound  over  the  barrier  and  reflections  from  the  other  barrier 
limit  the  reduction  of  the  energy  mean  level  at  this  receiver  to 
about  8 dB ( A) . Receiver  No.  2 is  only  partially  shielded  by  the 
reflective  barrier  against  noise  from  the  major  roadway  so  that 
the  noise  reduction  is  only  about  3 dB(A). 
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20 
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0. 
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10 

00. 
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20 

u 

1 
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300  . 

-5 
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-5 

5 

2 

550. 
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15 
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500. 

10 

6 

TEST 

CASE  2 

3 

2 
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100. 

30 
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100. 

20 

1 200. 

- 100. 

10 
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250  . 

40 
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250. 

40 

700. 
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40 
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300. 

40 

400. 

250. 
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RECEIVES  HEIGHT  ADJUSTMENT 
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STANDARD  DEVIATION  FOR  TRUCKS 
STANDARD  DEVIATION  FOR  CARS 
HEIGHT  ADJUSTMENT  FOR  CARS 
HEIGHT  ADJUSTMENT  FOR  TRUCKS 


L 


L 

L 

L 

240. 

T 

45.  2, 
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A 


R 
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TRAFFIC  NOISE  PREDICTION 


TEST  CASE  1 


program  initialization  PARAMETERS 

0 . 00000E+00 

1 RECEIVER  HEIGHT  ADJUSTMENT 

0 . 90000E+0 1 

2 number 

OF  FREQUENCY  BANDS 

0 . 25000E+0 1 

3 STANDARD  DEVIATION 

FOR 

TRUCKS 

0 . 35000E  + 0 1 

4 STANDARD  DEVIATION 

FOR 

CARS 

0 . 00000E+00 

5 HEIGHT 

ADJUSTMENT 

FOR 

CARS 

0 . 80000E+0  1 

6 HEIGHT 

ADJUSTMENT 

FOR 

TRUCKS 

ROADWAY  # 1 

number  OF 

VEH/H 

MPH 

TYPE  1 VEH 

1 

0. 1000E+04 

0 , 6000E+02 

NUMBER  OF 

VEH/H 

MPH 

TYPE  2 VEH 

1 

0. 1 000E+03 

0.5000E+02 

SOURCE  COORD 

IN  FT 

NUMBER  X 

Y 

Z 

1 0.1 000E+03 

0 , 2000E+03 

0.2000E+02 

2 0. 1 100E+04 

0 . 2000E+03 

0 , 0 0 0 0 E + 0 0 

3 0 . 1 6 00  E + 04 

0 . 0000E+00 

0 , 0000E+00 

ROADWAY  # 2 

NUMBER  OF 

VEH/H 

MPH 

TYPE  1 VEH 

1 

0 . 3000E  + 03 

0.4000E+02 

SOURCE  COORD 

IN  FT 

NUMBER  X 

Y 

Z 

1 0 . 7000E+03 

0 , 8000E+03 

0,  1 000E  + 02 

2 0 . 1 000E+03 

0.2000E+03 

0.2000E+02 

ABSORBING  STRIP  # 1 

(T) 

PT  # X 

Y 

Z 

WIDTH 

1 0.9000E+03 

0 , 3000E+03 

-0 , 5000E  + 0 1 

0 

. 2400E+03 

2 0.1  1 00E+04 

0 , 6000E+03 

-0 , 5000E+0 1 

RECEIVER 

RECEIVER  COORD  IN  FT 

NUMBER  X 

Y 

Z 

LBACK  SIDE  VLB 

1 0 . 5500E+03 

-0,  1 000E  + 03 

0,  1500E  + 02 

45.0  2.0 

2 0.5500E+03 

0 , 5000E+03 

0, 1000E+02 

50.0  3.0 

TEST  CASE  1 


RECEIVER  # 

XRC 

YRC 

ZRC 

1 

550,0 

-100,0 

15,0 

63 

40.7 

125 

53,0 

250  500 

57.9  62.9 

1000  2000 
62.5  60.3 

4000 

52,8 

8000 

39.4 

131 

LE(A) 

LNP 

L90 

L50 

L 1 0 

67.7 

78.4 

60.3 

65.7 

71,0 

INCH  . 

22.2 

27,8 

17,9 

20.7 

23,5 

2 

550.0 

500 

.0 

10,0 

63 

125 

250 

500 

1000  2000 

4000 

8000 

41.2 

53.2 

58.1 

63,1 

62,9  60,7 

53,2 

40.6 

LE  ( A ) 

LNP 

L90 

L50 

L 10 

68.0 

78,1 

61,1 

66.2 

71,3 

INCH. 

17.0 

19.4 

15.0 

16.2 

17,4 
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TRAFFIC  NOISE  PREDICTION 


TEST  CASE  2 
ROADWAY  # 1 


NUMBER 

OF 

VEH/H 

MPH 

TYPE  1 

VEH 

1 

0. 1000E+04 

0 , 6 000E  + 02 

NUMBER 

OF 

VEH/H 

MPH 

TYPE  2 

VEH 

1 

0. 1 000E+03 

0 , 5000E  + 02 

SOURCE  COORD 

IN  FT 

NUMBER 

X 

Y 

Z 

1 

0.  1 000E  + 03 

0 . 2000E  + 03 

0 , 2000E+02 

2 

0 . 1 1 00E  + 04 

0.2000E+03 

0 , 0000E+00 

3 

0 . 1 600E+04 

0.0000E+00 

0 , 0000E+00 

ROADWAY 

a 2 

NUMBER 

OF 

VEH/H 

MPH 

TYPE  1 

VEH 

1 

0 , 3000E+03 

0.4000E+02 

SOURCE  COORD 

IN  FT 

NUMBER 

X 

Y 

Z 

1 

0 . 7000E  + 03 

0 , 8000E  + 03 

0,  1000E  + 02 

2 

0. 1 000E  + 03 

0 , 2000E+03 

0 . 2000E+02 

BARRIER 

# 1 (A) 

barrier 

COORD  IN  FT 

NUMBER 

X 

Y 

Z 

1 

0 . 3000E+03 

0. 1 000E+03 

0 , 3000E+02 

2 

0 . 9000E+03 

0.  1 000E  + 0 3 

0 , 2000E+02 

3 

0. 1 200E  + 04 

-0.  1 000E  + 0 3 

0, 1 000E+02 

BARRIER 

# 2 (R) 

BARRIER 

COORD  IN  FT 

NUMBER 

X 

Y 

Z 

1 

0 . 400-0E  + 03 

0.2500E+03 

0 , 4000E+02 

2 

0 . 7000E  + 03 

0.2500E+03 

0 , 4000E+02 

3 

0.7000E+03 

0.3000E+03 

0 , 4000E+02 

4 

0 . 4000E+03 

0 , 3000E+03 

0.4000E+02 

5 

0 . 4000E+03 

0.2500E+03 

0 , 4000E+02 

ABSORBING  STRIP  # 1 

(T) 

PT  # 

X 

Y 

Z 

WIDTH 

1 

0 . 9000E+03 

0.3000E+03 

-0.5000E+01 

0 , 2400E  + 03 

2 

0. 1 1 00E  + 04 

0 . 6000E+03 

-0.5000E+01 

RECEIVER 

RECEIVER  COORD  IN  FT 

NUMBER 

X 

Y 

Z 

LBACK 

1 

0.5500E+03 

-0. 1000E+03 

0,  1500E  + 02 

45.0 

2 

0 . 5500E+03 

0.5000E+03 

0, 1 000E  + 02 

50.0 

STDEVLB 

2. 

3. 
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rj 


TEST  CASE  2 


RECEIVER 

# 

XRC 

YRC 

ZRC 

1 

550,0 

- 100 

.0 

15,0 

63 

36.  1 

125 

47,9 

250 

51.8 

500 

55.3 

1000 

53,0 

2000 

48,8 

4000 

39.7 

8000 

22,6 

INCR  . 

LE  ( A ) 
59.4 
14.0 

LNP 

75,6 

25,0 

L90 

46,8 

4.3 

L50 

54.8 

9.8 

L 1 0 

62,9 

15,4 

2 

553,0 

500 

.0 

10,0 

63 

38.1 

125 

49,5 

250 

54.1 

500 

59.1 

1000 

59.2 

2000 

57,0 

4000 

49,3 

8000 

36,9 

INCR  . 

LE  ( A ) 
64.4 
13.3 

LNP 
78,9 
20.  1 

L90 

53,4 

7.3 

L50 

60.7 

10.7 

L 1 0 

67,9 

14,1 
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Appendix  A: 


Memorandum  dated  2k  May  1972:  "Formula  for  Atmospheric 

Absorption  in  Computer  Program  for  the  Prediction  of  Road 
Traffic  Noise" 
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MEMORANDUM 


To:  Capt.  Wesler,  TSC  Cambridge 

From:  U.J.  Kurze 

Subject:  Formula  for  Atmospheric  Absorption  in  Computer  Program 

for  the  Prediction  of  Road  Traffic  Noise 

Date:  24  May  1972 


As  outlined  in  the  introduction  to  the  User's  Manual,  the 
computer  program  is  designed  to  predict  the  noise  mainly  within 
a distance  of  1000  ft  from  a highway.  At  distances  up  to  3000 
ft,  atmospheric  attenuation  plays  a minor  role  and  is  suffi- 
ciently taken  into  account  by  using  Eq . 42. 

Recent  attempts  to  predict  the  noise  in  rural  areas  at  dis- 
tances of  3 to  4 miles  from  a highway  yielded  results  for  the 
sound  pressure  level  which  look  too  high.  This  is  due  to  under- 
estimating the  atmospheric  absorption  in  the  500  Hz  octave  band 
which  is  the  most  typical  frequency  band  for  A-weighted  road 
traffic  noise.  In  the  frequency  band  below  1000  Hz,  relaxation 
processes  in  nitrogen  molecules  enhance  the  atmospheric  absorp- 
tion. 

This  enhancement  in  attenuation  at  lower  frequencies  as 
well  as  the  dominant  high  frequency  absorption  due  to  relaxation 
processes  in  oxygen  molecules  are  best  approximated  by  the 
formula 

Da  = 5.4  10"4  (2.35)n"5  r ' 

for  a temperature  around  6 8°P  and  relative  humidities  between 
50$  and  70$.  The  formula  yields  results  consistent  with  data 
published  by  C.M.  Harris,  "Absorption  of  Sound  in  Air  vs  Humid- 
ity and  Temperature,"  J.  Aooust.  Soc.  Amev.  40  (1966),  pp . 14S— 
159. 
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MEMORANDUM  TO  CAPT.  WESLER 
24  May  1972 
Page  2 

It  is  recommended  to  use  this  formula  in  the  computer  pro- 
gram and  to  replace  the  fourth  line  following  line  232  in  Sub- 
routine GEOMRY  by 

A = 10.**  (5.4  E-5  * 2.35  **  (IP-5)*  DN1 ) . 

Although  the  new  formula  is  aimed  at  improving  the  prediction  of 
traffic  noise  at  large  distances  from  a highway,  it  is  not  recom- 
mended to  use  the  computer  program  in  these  cases  without  careful 
considerations  of  possible  sources  of  large  errors,  such  as 
micrometeorological  conditions  and  ground  absorption. 
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